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A  unified  model  for  reverberation  and  submerged  target  scattering  in  a  stratified  medium  is 
developed  from  wave  theory.  The  advantage  of  the  unified  approach  is  that  it  enables  quantitative 
predictions  to  be  made  of  the  target-echo-to-reverberation  ratio  in  an  ocean  waveguide.  Analytic 
expressions  are  derived  for  both  deterministic  and  stochastic  scattering  from  the  seafloor  and 
subseafloor.  Asymptotic  techniques  are  used  to  derive  expressions  for  the  scattering  of  broadband 
waveforms  from  distant  objects  or  surfaces.  Expressions  are  then  obtained  for  the  scattered  field 
after  beamforming  with  a  horizontal  line  array.  The  model  is  applied  to  problems  of  active  detection 
in  shallow  water.  Sample  calculations  for  narrow-band  signals  indicate  that  the  detection  of 
submerged  target  echoes  above  diffuse  seafloor  reverberation  is  highly  dependent  upon  water 
column  and  sediment  stratification  as  well  as  array  aperture,  source,  receiver,  and  target  locations, 
in  addition  to  the  scattering  properties  of  the  target  and  seafloor.  The  model  is  also  applied  to 
determine  the  conditions  necessary  for  echo  returns  from  discrete  geomorphologic  features  of  the 
seafloor  and  subseafloor  to  stand  prominently  above  diffuse  seafloor  reverberation.  This  has  great 
relevance  to  the  geologic  clutter  problem  encountered  by  active  sonar  systems  operating  in  shallow 
water,  as  well  as  to  the  remote  sensing  of  underwater  geomorphology.  ©  2001  Acoustical  Society 
of  America.  [DOI:  10.1121/1.1339826] 

PACS  numbers:  43.30.Gv,  43.30.Hw,  43.30. Vh  [DLB] 


I.  INTRODUCTION 

A  common  problem  in  the  active  detection  and  localiza¬ 
tion  of  a  radar  or  sonar  target  arises  when  scattered  returns 
from  the  target  become  indistinguishable  from  returns  from 
randomly  rough  boundaries,  volume  inhomogeneities,  or  de¬ 
terministic  features  of  the  environment.  The  goal  of  the 
present  article  is  to  investigate  the  extent  to  which  environ¬ 
mental  reverberation  limits  the  ability  to  detect  and  localize  a 
target  submerged  in  an  ocean  waveguide,  where  methods 
developed  for  the  radar  half-space  problem  are  inapplicable 
due  to  the  added  complications  of  multi-modal  propagation 
and  dispersion. 

To  this  end,  a  unified  model  for  3-D  reverberation  and 
submerged  target  scattering  in  a  stratified  medium  is  devel¬ 
oped  from  wave  theory.  The  model  is  fully  bistatic  and  stems 
directly  from  Green’s  theorem,  since  it  generalizes  Ingenito’s 
approach1,2  for  harmonic  scattering  in  a  stratified  medium  by 
incorporating  stochastic  scatterers  and  time-dependent 
sources.  While  it  is  consistent  with  certain  narrow-band  re¬ 
sults  of  previous  “heuristic”3  derivations3-7  for  shallow  wa¬ 
ter  reverberation  measured  with  an  omni-directional  receiver 
that  are  based  on  the  work  of  Bucker  and  Morris,4  it  offers 
more  insight  and  generality  since  it  is  developed  from  first 
principles  with  explicitly  stated  assumptions.  For  example,  it 
clearly  obeys  reciprocity  for  source-receiver  locations 
within  a  layered  media,  which  is  important  in  properly  mod¬ 
eling  the  absolute  level  of  returns  from  targets  or  surfaces 
within  the  seafloor,  and  it  allows  absolute  comparison  be¬ 
tween  reverberation  and  deterministic  target  returns.  Such 
comparison  led  to  inconsistencies  in  previous  formulations 
as  noted  in  Ref.  3.  It  also  provides  analytic  expressions  for 
the  three-dimensional  (3-D)  field  scattered  bistatically  by 


both  stochastic  and  deterministic  objects  from  a  source  with 
arbitrary  time  function,  as  well  as  the  associated  spatial  and 
temporal  covariances.  This  enables  realistic  modeling  of  the 
moments  of  the  raw  reverberant  field  received  over  extended 
spatial  and  temporal  apertures  as  well  as  the  output  after 
subsequent  processing  with  standard  beamforming  and 
broadband  signal  processing  techniques.  In  the  present  ar¬ 
ticle,  applications  of  the  theory  are  restricted  to  systems 
which  employ  the  beamforming  and  temporally  incoherent 
processing  widely  used  in  narrow-band  signal  reception. 
Analytic  expressions  for  the  statistical  moments  of  the  scat¬ 
tered  field  are  obtained  directly,  but  can  also  be  obtained  by 
sample  averaging  over  realizations  by  Monte  Carlo  simula¬ 
tions,  as  for  example  is  done  for  rough  surface  scattering  in 
Ref.  8.  The  relative  merit  of  either  approach  depends  on  the 
relative  difficulty  in  evaluating  the  analytically  obtained  mo¬ 
ments  or  performing  the  Monte  Carlo  simulations  for  the 
given  problem.  The  analytic  approach  has  proven  to  be  more 
advantageous  and  insightful  for  the  illustrative  examples  of 
the  present  article. 

The  primary  motivation  for  developing  the  unified 
model  is  to  compare  the  absolute  level  of  target  echo  returns 
with  those  from  the  seafloor  and  to  investigate  how  these 
vary  in  both  absolute  and  relative  level  as  a  function  of  water 
column  and  sediment  stratification,  receiving  array  aperture, 
and  source,  receiver,  and  target  locations  in  a  shallow  water 
waveguide.  Another  major  focus  of  the  present  article  is  to 
investigate  the  manner  in  which  scattering  from  both  ex¬ 
tended  geomorphologic  features  and  randomly  rough  patches 
of  the  seafloor  and  subseafloor  contribute  to  measured  rever¬ 
beration.  The  latter  typically  makes  up  the  diffuse  reverber¬ 
ant  background,  which  has  an  expected  intensity  that  decays 
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in  time  after  the  arrival  of  the  direct  signal  waveform.  The 
former  typically  leads  to  geological  clutter,  which  is  defined 
as  any  set  of  acoustic  returns  from  the  seabed  that  stand 
significantly  above  the  diffuse  and  temporally  decaying  re¬ 
verberation  background.  Geological  clutter  is  a  primary 
problem  in  active  sonar  operations  in  shallow  water.  This  is 
because  the  clutter  can  be  confused  with  or  camouflage  re¬ 
turns  intended  from  a  submerged  target. 

A  goal  of  this  work  is  then  to  determine  plausible  physi¬ 
cal  mechanisms  for  geological  clutter  by  use  of  the  unified 
model.  To  do  so,  both  first-order  perturbation  theory  and 
empirical  Lambert-Mackenzie9  models  are  used  to  describe 
stochastic  scattering  from  a  randomly  rough  seafloor.  These 
models,  together  with  deterministic  models  for  scattering 
from  seafloor  features,  are  used  to  determine  scenarios  where 
geologic  clutter  is  significant.  The  present  focus  is  on  sea¬ 
floor  features  with  mean  surfaces  that  are  finite  and  inclined 
with  the  reflection  properties  of  the  layer  to  which  they  be¬ 
long,  since  such  features  are  ubiquitous  in  continental  shelf 
waters.  While  a  large  literature  exists  for  scattering  from  2-D 
features  in  a  waveguide,  the  focus  of  the  present  work  is  on 
scattering  from  3-D  features  in  a  waveguide.  Apparently,  the 
only  previous  work  on  deterministic  scattering  from  3-D  sea¬ 
floor  features  in  a  waveguide  has  been  for  acoustically  com¬ 
pact  (ka<g  1)  proturbances  on  perfectly  reflecting  bottoms,10 
but  compact  targets  are  too  weak  to  comprise  geological 
clutter  in  a  long-range  active  sonar  system  and  are  not  rel¬ 
evant  to  the  present  analysis.  A  review  of  the  general  litera¬ 
ture  on  3-D  scattering  in  an  ocean  waveguide  is  given  in  Ref. 
11.  All  illustrative  examples  in  the  present  article  employ 
time-windowed  cw  source  waveforms  and  monostatic  geom¬ 
etries  to  investigate  the  central  detection  issues  with  as 
simple  an  approach  as  possible. 


II.  THE  UNIFIED  MODEL 

A.  The  single-scatter  approximation  for  3-D  scattering 
from  an  object  of  arbitrary  shape  in  a  layered 
medium 

A  number  of  simplifying  conditions11  that  apply  to  a 
wide  variety  of  active  sonar  problems  in  the  ocean  enable  the 
field  scattered  from  an  object  submerged  in  a  stratified  me¬ 
dium  to  be  approximated,  from  Green’s  exact  theorem,  as  a 
linear  function  of  the  object’s  plane  wave  scatter  function.12 
The  plane  wave  scatter  function  S(0,cf>\  0, ,  (/>,)  at  frequency/ 
is  defined  in  Appendix  A,  where  its  relationship  to  Green’s 
theorem  and  the  traditional  target  strength  and  scattering 
strength  measures  of  ocean  acoustics  is  explained. 

To  formulate  the  unified  model,  it  is  convenient  to  ini¬ 
tially  follow  Refs.  1,  2,  and  11  by  placing  the  object  centroid 
at  the  center  of  all  coordinate  systems.  The  source  coordi¬ 
nates  are  then  defined  by  (x0,yQ,Zo)>  receiver  coordinates  by 
(x,y,z),  and  coordinates  on  the  surface  of  the  object  by 
( xt,yt,zt )  where  the  positive  z  axis  points  downward  and 
normal  to  the  interfaces  between  horizontal  strata.  Spatial 
cylindrical  ( p,6,z )  and  spherical  systems  (/-,  0,  </>)  are  de¬ 
fined  by  x  =  r  sin  #cos  </>,  y  =  r  sin  6  sin  </>,  z  =  rcos9,  and 
p2=x2+y2.  The  horizontal  and  vertical  wave  number  com¬ 
ponents  for  the  «th  mode  are  respectively  £„  =  k  sin  a„  and 


y„  =  k  cos  an ,  where  k2  =  £2  +  y2t  and  the  wave  number  mag¬ 
nitude  k  equals  the  angular  frequency  oj  divided  by  the  sound 
speed  c  in  the  target  layer. 

The  spectral  component  of  the  scattered  field  for  a  time- 
harmonic  source  of  frequency  /  at  r0  and  a  receiver  at  r  then 
becomes 


oo  oo 

®,(r|r0)~  2  2  <J>‘m'n)(r|r0),  (la) 

m= 1 n= 1 


where 


<£*m’")(rl r0)  =  ~Y  [^m(rMn(ro)5'( 77-  a,„ ,  a„  ,(f>Q+i r) 

~  B  ( r ) A  „  ( i r0 ) , S  ( am  ,  / ;  an  ,cf>0  +  tt) 
-Am(r)R„(r0)S( n t- am  ,  <f>\ tt- a„  ,  <f>0+  tt) 
+  Bm{r)Bn(r0)S(am ,  tt-  an  ,cf>0+i r)], 

(lb) 

and 


AJr)  — (8^mp)  mum{z)N mel(^p+y«<D  17/4),  (2a) 


Bm(r)  =  d{0j (2b) 


A„(ro)=^(877^p0)“1/2M„(z°)A„-ei(^0+^-ir/4), 


(2c) 


fi,Xro)=^(8^„p0)-1/2«,/z0)<e'(^0-2'"D-W4), 

(2d) 

are  the  down-  and  up-going  plane  wave  amplitudes  in  the 
layer  of  the  object,  D  is  the  depth  of  the  object  center  from 
the  sea  surface,  d(z)  is  the  density  at  depth  z,  and  ujz)  are 
the  mode  functions.  The  product  of  e  and  the  right- 

hand  side  of  Eq.  (la)  yields  the  time-harmonic  scattered 
field.  The  mode  functions  are  normalized13  according  to 

f”  <{z)un{z) 

Snm  J-D  d(z)  dz’  (3) 

and  must  be  decomposable  into  up-  and  down-going  plane 
waves  via 

un(z)  =  N~eiy"U+D)-N^e-i^z+D'>  (4) 

in  the  layer  of  the  object,  where  N~  and  A4  are  the  ampli¬ 
tudes  of  down-  and  up-going  plane  waves  in  this  layer.  In  a 
Pekeris  waveguide,  for  example. 
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n„  =n : 


V2  1/  sin  2  y„H\ 
2  i  d\H  2  yn  ) 


1  •  9  U  ~  1/2 

1  sin  ynH 

-I -  —  (5) 

<4  1 

where  r/;,  and  ch  are  the  density  and  sound  speed  of  the 
bottom.  If  the  mode  functions  are  specified  at  any  two 
depths,  z.  1  and  z2  within  the  target  layer,  the  down-  and  up- 
going  plane  wave  amplitudes  can  be  readily  obtained  as 


n;  ' 

e‘Vn^l+D) 

e~iyn(z  i+D) 

-1 

«n(Zl) 

-N+ 

1  n 

eiy„(z2+D) 

g-'JniZl  +  D) 

un(z2). 

Equations  (l)-(6)  for  the  scattered  field  from  an  arbitrarily 
shaped  object  in  a  waveguide  differ  from  Ingenito’s 
formulation1  in  a  number  of  ways.  The  most  substantial  dif¬ 
ference  is  that,  by  inclusion  of  Eqs.  (4)  and  (6),  they  explic¬ 
itly  show  how  the  scattered  field  for  an  arbitrarily  shaped 
object  can  be  computed  in  a  stratified  medium.  Ingenito  also 
defines  the  plane  wave  scatter  function  differently  than  most 
standard  texts  by  describing  the  incident  plane  wave  in  terms 
of  the  direction  it  comes  from  rather  than  the  direction  it 
goes  to.  The  latter,  standard  approach,  is  adopted  here. 

Finally,  the  more  standard  mode  function  normalization 
of  Ref.  13  is  adopted  here,  so  that  Eq.  (1)  obeys  reciprocity 
as  defined  in  Appendix  A2  of  Ref.  13,  so  that 
^(zo)9T5(r|r0)==£/(z)<I>s(r0|r).  Satisfaction  of  reciprocity 
becomes  important  for  an  approach  if  it  is  to  yield  accurate 
estimates  of  the  scattered  field  when  the  source,  receiver,  and 
target  are  in  layers  that  have  significantly  different  densities, 
and  is  a  natural  consequence  of  the  use  of  Green’s  theorem 
in  the  present  formulation,  but  has  been  left  unaddressed  as 
an  issue  in  previous  heuristic  reverberation  formulations. 3-7 
The  issue  becomes  of  practical  concern  in  modeling  the  level 
of  returns  from  targets  or  surfaces  buried  in  the  seafloor  from 
sonar  systems  operating  in  the  water  column  above. 

A  more  general  expression  than  Eqs.  (1)— (6),  for  the 
scattered  field  from  an  arbitrarily  shaped  object  in  a  stratified 
medium,  is  given  in  Refs.  2  and  1 1  in  terms  of  wave  number 
integrals.  As  noted  in  Ref.  11,  the  more  general  wave  num¬ 


ber  formulation  is  valid  when  (1)  the  propagation  medium  is 
horizontally  stratified  and  range  independent;  (2)  the  object 
is  contained  within  an  iso-velocity  layer;  (3)  multiple  reflec¬ 
tions  between  the  object  and  waveguide  boundaries  make  a 
negligible  contribution  at  the  receiver;  and  (4)  the  range  from 
the  object  to  source  or  receiver  is  large  enough  that  the  scat¬ 
tered  field  can  be  expressed  as  a  linear  function  of  the  ob¬ 
ject’s  plane  wave  scatter  function.  All  of  these  conditions 
then  must  be  satisfied  for  Eqs.  (1)— (6)  to  be  valid,  with  one 
additional  constraint.  The  ranges  involved  must  be  large 
enough  that  the  Green’s  functions,  from  source-to-target  and 
target-to-receiver,  are  accurately  approximated  as  sums  of 
discrete  modes.  The  latter  is  an  expected  consequence  of  the 
Riemann-Lebesgue  lemma. 14  The  present  formulation  and 
its  spectral  equivalent  have  been  implemented  for  target  scat¬ 
tering  in  a  waveguide  over  the  full  360-degree  span  of  bi¬ 
static  angles  in  Refs.  1,  2  and  11.  It  is  noteworthy  that  this 
formulation  includes  the  scattering  of  evanescent  waves  by 
analytic  continuation  of  the  scatter  function,  as  has  been  pre¬ 
viously  discussed  and  implemented  in  Refs.  2  and  1 1,  as  well 
as  in  Ref.  8  which  uses  a  formulation  similar  to  Ingenito’s. 


B.  The  field  scattered  from  general  stochastic  targets 

By  allowing  the  scatter  function  for  the  object  to  be  a 
random  variable,  the  single-scatter  formulation  of  the  previ¬ 
ous  section  applies  to  the  more  general  problem  of  scattering 
from  a  stochastic  target  submerged  in  a  waveguide.  This  ap¬ 
proach  is  particularly  valuable  in  modeling  scattering  from 
targets  of  unknown  shape  or  orientation,  randomly  rough 
surface  interfaces,  or  stochastic  volume  heterogeneities,  all 
of  which  can  contribute  significantly  to  the  reverberant  field 
measured  in  shallow  water. 

The  moments  of  the  scattered  field  can  be  derived  ana¬ 
lytically  to  determine  its  expected  behavior.  The  mean  field, 
for  example,  becomes 

00  00 

<3>,(r|r0)>~  S  2  <<^m’n)(r|r0)>,  (7) 

m= 1 n = 1 

where 


^(m,«)(r|r0))=  —[Am(r)An(r0){S(TT-am,(f>\an,<f>0+Tr))-Bm{r)An(r0){S(am,<f>\an,<f>0+TT)) 

-Am(r)Bn(r0)(S (it- am ,  <£;  77-  a„ ,  (f>Q+  77))  +  B,„(r)B„(r0)(S(  am ,  <f>;  77-  a„  ,  <^0+  77))] 
while  the  mutual  intensity  of  the  field  scattered  for  receivers  at  r  and  r'  becomes 

co  00  CO  co 

<O>.s(r|r0)D?(r'|r0))~2  S  S  S  <<^n)(r]ro)3>(m'’"9  Vlr0)>, 

m=l  n=  1  m'  =  1  „f  =  1 

where 


(8) 

(9) 


[ Am(r)An(ro)A*,(r')A*(ro)(S('n—am,0;an,</>o+TT)S*(TT-am,,cf>';an,,</>o+'n j) 

-Am(r)An(r0)B*,(r')A*  (r0)(S(n-  a,  „,</>;<?„,  c/>0+ ir)S*(am,  ,(/>'  </>„+  77)) 

— 4m(r)A„(r0)A*,(r')5*  (ro)<5(77—  a, „,4>\an,  f0+ir)S*(n- am,  ;  17 77)) 
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+  Am(r)An(r0)B*,(r')B*  (r0)(S(-n—a,  n,ff,an,  <f>0+  Tr)S*(am,  ,(/>'  -,tt-  an,  ,c/>0+  tt)) 

~ B m(r)An(r0)A*,(r' )A*,(r0)(S(am an ,(/>0+  i am,  ,</>' \an,  ,(f>0+  tt)) 

+  Bm(r)An(r0)B*,(r')A*,  (r0)(S(a,  „,</>;<*„,  4>0+ ir)S*(am, ,  <t>o+ n)) 

+Bm(r)An(r0)A*,(r')B*(r0)(S(am,<l>-,an,<f>0+Tr)S*(Tr-am,,cl>'-,'n—an,,<l>0+Tr)) 
-Bm(r)An(r0)B*,(r')B*,  (r0)(S(a,  4>0+  ir)S*(am,  ,tf>'  ;tt-  an,  ,(f>0+  tt)) 

-Am(r)Bn(r0)A*,(r')A*(r0)(S(ir-a!m,</>-,TT- an,(/>0+  am,  ,c/>'  ;an,  ,(/>0+  tt)) 

+Am(r)Bn(r0)B*,(r')A*(r0)(S('TT-am,(f>;'n—  an,(f)0+  n)S*(am, ,(/}' \an,  ,<f>0+  it)) 
+Am(r)B„(ro)A*,(r')5*  (r0)(5(T7-  am,  </>;  it- an,(/>0+  n)S*(n~  am,  ,</>'■, tt- an,  ,<f>0+  tt)) 
-Am(r)B„(r0)B*,(r')B*, (r0)(S(ir- am a„  ,</>„  +  ir)S*(am,  ,</>'; tt- an,  ,</>0  +  tt)) 

+  Bm(r)B„(r0)A*, (r')A*,(r0)(S(a:m, <t>;Tr-  a„,<f>0+  ir)S*(Tr- am,  ,4>0+  tt)) 

-Bm(r)B„(r0)B*,(r')A*,  (r0)(S(  am,ff tt-  a„,<f>0+  ir)S*(am,  ,(/>' -,an,  ,cf>0+  tt)) 
-Bm(r)Bn(r0)A*,(r')B*,(r0)(S(a:m,<t>;Tr-all,<f>0+Tr)S*{Tr-am,,(l>'-,Tr-an,,<t>0+Tr)) 

+  Bm(r)Bn(r0)B*,(r')B*{T0)(S(am,</>;ir-an,<ft0+-ir)S*(am,,4>'\ir-an,,<f>0+ir))].  (10) 


The  spatial  covariance,  or  cross  spectral  density,  of  the  scat¬ 
tered  field, 

($J(r|ro)$f(r'|ro))-(O)J(r|ro))(<Df(r']r0)) 

00  00  00  00 

H  ($«"'")(r|r0)®("'-"')*(r'|r0)) 

m=\  n=  1  m'  =  1  „'=1 

-(c&^(r|r0))(<Dl'"'-"'^(r'|r0)),  (11) 

follows  directly  from  Eqs.  (7)-(10),  and  has  been  imple¬ 
mented  for  fluctuating  targets  submerged  in  an  ocean 
waveguide.15 

C.  When  scattering  statistically  decorrelates  the 
waveguide  modes 

Many  useful  scattering  properties  of  a  random  target  are 
described  by  the  covariance  of  its  scatter  function,  which 
couples  incident  modes  with  n  and  n '  subscripts  together 
with  scattered  modes  with  m  and  m '  subscripts  in  the  cross 
spectral  density  of  Eq.  (11).  Under  many  situations  of  prac¬ 
tical  interest  in  the  scattering  from  fluctuating  targets  or  ran¬ 
domly  rough  surfaces,  the  n  and  n '  incident  modes  decouple 
as  do  the  scattered  m  and  m '  modes. 

Assume  that  the  random  target’s  characteristic  dimen¬ 
sion  L  is  large  compared  to  the  wavelength  \  and  that  for  any 
realization  it  has  complicated  structure  with  significant  varia¬ 
tions  on  the  wavelength  scale.  These  assumptions  are  used 
extensively  in  radar  to  describe  fluctuating  targets,16’17  such 
as  aircraft  of  unknown  shape  and  orientation,  and  in  radar, 
statistical  optics,  and  acoustics  to  describe  scattering  from 
randomly  rough  surfaces.18  20  They  typically  lead  to  circular 
complex  Gaussian  random  (CCGR)  fluctuations  in  the  scat¬ 
tered  field  over  different  statistical  realizations  of  the 


target. 16-21  When  the  target  is  a  randomly  rough  surface, 
these  assumptions  are  equivalent  to  restricting  the  correlation 
length  of  the  surface  roughness  to  be  much  smaller  than  the 
dimension  L  of  the  targeted  surface.  In  this  case  the  scattered 
field  is  usually  termed  diffuse  or  incoherent  “because  of  its 
wide  angular  spread  and  lack  of  phase  relationship  with  the 
incident  wave.”18 

For  any  realization  of  such  a  random  surface  or  target, 
say  for  example  the  £th  realization,  its  scatter  function 
will  be  highly  oscillatory  in  magnitude  and  phase 
over  both  incident  O and  scattered  {l  =  (6,<j>)  di¬ 
rections  and  will  have  a  complicated  lobe  pattern.  This  lobe 
pattern  will  vary  significantly  over  random  realizations  of  the 
surface  or  target  due  to  changes  in  constructive  and  destruc¬ 
tive  interference,  making  the  scatter  function  representing  all 
realizations  5(fi ;£!,-)  a  random  variable  with  zero  expected 
value,  (5(0; ft,-))  =  0.  Due  to  the  complicated  structure  of 
the  target,  the  angular  width  of  any  lobe  in  the  scatter  func¬ 
tion  will  be  on  the  order  of  the  minimum  width  set  by  dif¬ 
fraction  of  X/L,  which  is  small  by  assumption.  The  second 
moment  of  the  scatter  function  (5(n;n,)5*(n';n,-))  for 
incident  angles  il, ,  il-  and  scattered  angles  ft,  ft'  can  be 
thought  of  as  an  ensemble  average  of  Sk{fl',fli)S*(fl' ’,0,-) 
over  realizations  of  the  random  surface.  This  product  will 
oscillate  about  zero  across  realizations  k  for  0,  —  ft' 
>  k/L  and  |  ft  —  ft '  \  >  k/L  so  that  the  ensemble  average  will 
tend  to  zero  due  to  term  by  term  cancellation.  Here  k/L 
behaves  as  an  angular  correlation  width  of  the  scatter  func¬ 
tion  over  variations  in  both  incident  and  scattered  angle.  The 
product  will  tend  to  become  positive  definite  for  it  ,■  —  ft  ■ 
=£ k/L ,  | ft  —  n'|=£\/L,  however,  so  that  the  ensemble  aver¬ 
age  will  approach  (|5(ft;ft,)|2). 

Applying  this  reasoning  to  stochastic  targets  or  surface 
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patches  in  a  waveguide,  the  covariance  of  the  scatter  function  can  be  written  as 

(S(  am  ,  (f>\ a„  ,  (f>0  +  v)S*  ( am  a„ ,,  (f>0  +  tt))  -  (S(  am  ,  0;  a„  ,  cf>0  +  it)  )(S  *(  am  a„ ,,  (f>0  +  tt)) 

=  8nn’Smm'Cmn{  r,r'|r0),  (12a) 

where 

Cm„(r,r'\rQ)  =  {S(am,(/)-,an,(f>0+TT)S*(am,(f>';an,(t>0+'Tr))-(S(an,,cl>-,an,(f>0+Tr))(S*{am,(f>’-,an,(t)0+TT)),  (12b) 


when  the  angular  separation  |  an  —  an ,  between  modes  n  and 
n '  is  greater  than  the  angular  correlation  width  X/L  set  by 
diffraction.  In  practical  scenarios,  this  criterion  will  apply 
most  stringently  to  the  dominant  modes,  where  full  decorre¬ 
lation  will  occur  for  a„  ,  ,  —  an  on  the  order  or  greater  than 
X/L.  For  trapped  modes,  |  an ,  —  an  is  less  than  the  critical 
angle  of  the  bottom. 

Insertion  of  Eq.  (12a)  into  Eq.  (11)  then  leads  to  a  great 
simplification  in  the  spatial  covariance  of  the  scattered  field 

(<DJ(r|r0)<i>f(r'|r0))-(<DJ(r|r0))(d)f(r'|r0)) 

00  00 

-EE  (^(rMa^Vlro)) 

m  =  1  n  —  1 

—  <<F^n)(r|r0)><cF^V'|r0)),  (13) 

where  the  quadrupole  modal  sum  reduces  to  a  double  modal 
sum. 


D.  The  field  scattered  from  a  randomly  rough  or 
inhomogeneous  seabed 

When  the  scattering  is  due  to  a  randomly  rough  seafloor 
patch,  a  great  simplification  occurs  in  the  form  of  the  mutual 
intensity.  Of  the  16  parenthetical  terms  of  Eq.  (10),  only  1 
represents  a  down-going  incident  wave  coupling  to  an  up- 
going  scattered  wave.  Accordingly,  the  field  scattered  into 
the  waveguide  from  a  randomly  rough  seafloor  patch  of  area 
A  A,  as  defined  in  Appendix  A,  must  have  mutual  intensity 
given  by  Eq.  (9)  with 

((p('«.")(r|r0)cp<'",’n,)*(r'|r0)) 

/4tt\2 

=  (  ~r)  )A*'(ro) 

X{S(am,(f>;an,(f>Q+Tr)S*(am,  ,(f>0+ir))], 

(14) 

and  cross  spectral  density  given  by  Eq.  (11)  with 


^<p(m,n)(r|ro)<p(m'>"')*(r'|r0))  — (<&l'n’'4(r|r0))(<Fjm,’",l*(r'|r0)) 


4t t 


=  “H  fim(r)A„(r0)fi*,(r')A*  (r0){{S{am,(f>;an,(t>0+ Tr)S*{am,  ,(f>'\an,  .^q+tt)) 


-  ( S ( am  ,  0;  an  ,  (/>o  +  tt) ) ( S * ( a,„ , ,  4> ' ;  an , ,  <£0  +  tt ) ) }, . 


(15) 


The  cross  spectral  density  is  more  useful  in  describing  the 
stochastic  scattering  properties  of  a  randomly  rough  seafloor 
patch  than  the  mutual  intensity  because  deterministic  effects, 
such  as  specular  reflection,  coherent  beaming,  and  forward 
scattering,  are  removed  with  the  expected  field.  In  diffuse 
surface  scattering  problems,  where  the  surface  scattering 
patch  must  be  much  larger  than  the  wavelength,  the  expected 
value  of  the  scattered  field  is  typically  negligible  away  from 
the  specular  direction  due  to  random  interference.  The  cross 
spectral  density  and  mutual  intensity  then  become  effectively 
indistinguishable. 

With  the  assumption  of  diffuse  scattering  described 
in  the  previous  section,  which  is  supported  by  a  large 
amount  of  experimental  evidence,18,22  application  of  Eq. 
(12a)  yields 


(<D^,n)(r|r0)0)^',n')*(r'|r0)) 

-(cF<"1,")(r|r0))(c&(m',"'^(r'|r0)) 

=  Smm'Snn'  ,Cmn(r,r'\r0)Bm(r)An(r0)Bl,(r')A*, (r0) 


which  leads  to  a  great  simplification  in  the  scattered  field 
covariance: 

(<Ds(r|r0)<Ds*(r'|r0))-((Ds(r|r0))(a)f(r'|ro)) 

I  4tt\2  ”  " 

=  -H  E  E  Bffl(r)fi*(r')lA»(r0)|2Cm„(r,r'|r0). 

\  k  )  m= 1 n= 1 

(17) 
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If  the  seafloor  is  taken  as  an  aggregate  of  range-dependent 
scattering  patches  and  if  these  are  each  small  enough  to  have 
little  effect  on  the  mean  forward  field,  Eq.  (17)  provides  a 
good  approximation  to  the  cross  spectral  density  after  for¬ 
ward  propagating  through  this  mildly  range-dependent  wave¬ 
guide  and  then  scattering  in  a  specified  direction  from  the 
given  patch. 

It  is  noteworthy  that  a  formulation  in  terms  of  wave 
number  integrals  is  not  convenient  in  describing  the  statisti¬ 
cal  equipartition  of  energy  associated  with  diffuse  scattering 
since  the  modes  are  the  entities  that  describe  the  system’s 
degrees  of  freedom  rather  than  the  wave  number  compo¬ 
nents. 


E.  General  saddle  point  approximation  for  the 
scattered  field  in  time  from  a  distant  object 

For  a  source  with  general  time  dependence 

q(t)=  r  Q(f)e~i2vft df,  (18) 

J  —CO 

and  spectrum 

Q(f)=  r  q(t)ei2”"dt,  (19) 

J  —00 

the  scattered  field  as  a  function  of  time  from  an  object  with 
center  at  the  origin  becomes 

*,(rkolO=  P  Q (/) cI\s ( r | r0 ) e ~ i2lrft df,  (20) 

where  Eq.  (1)  can  be  rewritten  as 

<t\5(r|r0)=2  S  Am„(rlr0 (21) 

m  =  1  n  =  1 

so  that 

r  oo 

*,(r|r0|f)=2  2  2(/)Am«(r|r0 ,f)eip^mn^>  df , 

m=  1  n=\  J  —oo 

(22) 

where 

KJr\r0,f)  =  ^sm’n)(r\r0)e^P^‘P^,  (23) 

and 


Ltf) 


Zm-lTtf-- 


(24) 


By  application  of  the  saddle  point  method,  for  large  p,  as¬ 
suming  p0/p  and  f/p  are  fixed. 


^s(rlro]0-2  2  2 

/  =  1  m—  1  «=1 


277 


P<P"{flmn) 


dl/'mntf ) 


df 


=  0, 


f=fl„ 


(26) 


and  choosing  the  complex  roots  that  lead  to  a  finite  solution 
of  Eq.  (24)  as  p  and  p0  increase,  where  /  is  the  index  for 
multiple  roots  given  n  and  m.  Such  solutions  do  not  exist 
before  the  wave  packet  has  arrived  or  after  it  has  passed  the 
receiver.  Equation  (24)  can  also  be  written  more  conve¬ 
niently  as 


^s(r|r0|f)~2  2  2  Qiflmn ) 

1=  1  m=  1  n  =  1 


2tt 


Pf'fflmn) 


X  ( ” 'm  )(  r  I  r0)  I  e~  i2/K,flmn  +  (  ±  « 

s  J  J  Imn 

(27) 

where  the  frequencies  //m„  for  each  n  to  m  mode  conversion 
must  be  evaluated  at  each  time  and  source  and  receiver 
range.  Equation  (27)  is  also  obtained  if,  analogously,  p0  is 
made  large  and  r/p0  held  fixed  in  Eqs.  (22)-(24).  Typically, 
both  p  and  p0  will  be  sufficiently  large  for  the  saddle  point 
method  approximation  to  hold  whenever  the  modal  formula¬ 
tion  of  Sec.  II  A,  which  also  requires  large  p  and  p0 ,  is  valid. 

The  covariance  of  the  scattered  field  at  time  t  from  a 
distant  stochastic  target  then  becomes 

(^J(r|r0|f)^f(r'|r0|t))-(^i(r|r0|t))(^f(r'|r0|t)) 


-22  2  2  2  2  Q{flmn)Q*if Vm>n’) 

1=1  m  =  1  n=  1  P  =  1  m'  =  i  „'  =  1 


277 


X- 


V 1 PP'  f"(flmn)f"Xfl'm'n') 

Imn  —  f i  •  m  >  n  1 7r^)  Imn  —  '  77’/4)  / '  m'  n’ 


-(3>>"’(r|r0)|/=/w 


(28) 

by  substitution  of  Eq.  (la)  as  appropriate  for  the  given  scat¬ 
tering  scenario. 

It  can  now  be  seen  that  the  basic  equation  of  Ref.  3,  Eq. 
(9),  which  appears  without  stated  restrictions,  is  not  gener¬ 
ally  valid,  except  under  certain  narrow-band  conditions. 
Also,  the  present  analysis  indicates  that  the  group  velocity 
cannot  generally  be  treated  as  a  frequency-independent  quan¬ 
tity  as  it  is  in  the  development  of  Ref.  3,  where  a  number  of 
narrow-band  assumptions  have  apparently  been  made  im¬ 
plicitly,  as  may  be  seen  by  also  consulting  Refs.  7  and  23,  for 
example. 


F.  An  absolute  reference  frame 


xQifi  mn  )A  n  in  (r|r  0,/z  mn) 

Xe'Pt  ( flmn )  +  (  ±  I  ■ ir/4)lmn ;  (25  ) 

where  the  relevant  saddle  points  fimn  are  determined  by  solv¬ 
ing  the  equation 


To  compute  reverberation  from  wide  and  heterogeneous 
areas  of  seafloor  or  a  number  of  distributed  scatterers,  it  is 
convenient  to  recast  Eq.  (1)  in  terms  of  an  absolute,  rather 
than  target-centered,  spatial  coordinate  system.  Let  this  sys¬ 
tem  be  defined  by  coordinates  R  =(X,Y,Z)  whose  axes  are 
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parallel  to  those  of  the  x,y,z  target-centered  system,  where 
the  positive  Z  axis  is  again  downward  pointing,  but  whose 
origin  lies  at  the  ocean  surface,  for  example,  where  z=~D 
in  the  target-oriented  frame.  In  this  more  general  frame,  the 
source  position  is  defined  by  R0  =  (X0,y0,Z0),  the  receiver 
positions  by  Rr=(Xr,Yr,Zr)  and  R'  =  (X'r  ,Y'r  ,Z'),  and  the 
center  of  a  given  scattering  patch  by  R=  (X,  Y,Z),  where,  for 
example,  X  =  R  sin  i9xos  (p,  Y  =  R  sin  fl  sin  <p,  Z=R  cos  fl, 
and  R2  =  X2  +  Y2  +  Z2.  The  origin  of  all  these  coordinate  sys¬ 
tems  are  colocated  and  the  axes  are  parallel.  Spatial  coordi¬ 
nates  are  translated  from  the  target-oriented  to  the  absolute 
frame  by  substituting  r=R,.  — R,  r'=R'— R,  and  r0  =  R() 
-R  in  Eqs.  (1)  and  (2).  This  leads,  for  example,  to 


_ ium(Z,-Z)Nm _ 

d{Z)  ( 8  t r£„  V(. X-X)2+(Y-Y )2) 1/2 


Xei(^(Xr-x)2+(,Y-Y)2+ymZ-M)t  (29a) 


sm(R,-R) 


_ i“m(Zr-Z)N+ _ 

d(Z)(  8  7t£„  ^(X,-X)2+(Y-Y)2) 1/2 


Xei(Zm^(Xr-X)2+(Yr-Y)2~YmZ-TT/4)  (29b) 


reverberation  to  be  diffuse,  the  scattering  region  that  contrib¬ 
utes  to  the  intensity  measured  at  a  given  instant  must  be  large 
compared  to  the  mean  wavelength.  This  region,  referred  to  as 
the  system  resolution  footprint,  will  be  considerably  smaller 
for  data  beamformed  with  a  high-resolution  array  than  for 
data  received  by  an  omni-directional  receiver.  The  second 
component,  known  as  clutter,  is  here  defined  as  any  discrete 
temporal  event,  caused  by  an  anomalous  scatterer,  that 
stands  significantly  above  the  diffuse  reverberation  back¬ 
ground.  Here  “significantly  above”  means  much  more  than 
one  standard  deviation  in  sound  pressure  level.  For  certain 
systems  that  employ  high-resolution  temporal  processing  and 
operate  in  weakly  dispersive  waveguides,  there  may  be  no 
diffuse  component  to  the  reverberation.  In  this  case  coherent 
temporal  oscillations  may  be  found  in  reverberant  intensity 
measurements24  that  are  due  to  modal  interference  as  noted 
by  Ellis.3  Lepage  has  recently  investigated  similar  coherent 
effects  under  a  narrow-band  approximation  for  an  omni¬ 
directional  receiver.7 

When  the  single-scatter  approximation  is  valid,  the  total 
reverberant  field  measured  at  any  time  t  for  a  time-harmonic 
source  is  simply  the  sum  of  the  scattered  fields  from  all  en¬ 
vironmental  scatterers 


by  making  the  substitutions 

x0=(X0-X),  y0=(Y0-Y), 

x=(X-X),  y  =  (Y-Y), 

Po=\I(X0-X)2  +  (Y0-Y)2,  (30a) 

p=^X-X)2  +  (Y-Y)2, 


<Dr(Rr,R0)e-‘'2^ 

=  e~i2nf,J  J  J <t>s(Rr -R\R0-R)dXdYdZ.  (32) 

v 

For  a  source  with  general  time  dependence  q(t),  the  total 
reverberant  field  becomes 


—  Z0  Z,  z  —  Zr  Z, 
in  Eqs.  (1)  and  (2),  where 

cos  <£0  =  x0/p0,  sin  </>0  =  y0/po, 

,  ,  ■  ,  ,  (30b) 

cos  <p=x/p,  sin  <p=y/p. 

It  must  be  stressed  that  the  plane  wave  amplitudes  and 
vertical  wave  numbers  are  evaluated  in  the  layer  of  the  scat¬ 
tering  patch.  The  covariance  of  the  scattering  function  for  a 
given  patch  in  the  absolute,  rather  than  object-oriented  frame 
then  becomes 

cmn  ( R,  -  R,  R,'  -  R|  R()  -  R)  =  cmn{  r,  r '  |  r0 ) .  (31) 

III.  SHALLOW  WATER  REVERBERATION 
A.  Reverberation  in  time 

Reverberation,  as  measured  with  an  active  sonar  system, 
is  taken  to  be  any  and  all  echoes  returning  from  the  environ¬ 
ment  rather  than  the  intended  target.  The  characteristics  of 
reverberation  then  depend  not  only  on  the  environment  but 
also  the  geometry  of  the  source  and  receiver  as  well  as  the 
signal  waveform.  In  field  measurements,  reverberation  is 
measured  as  a  function  of  time.  It  can  often  be  decomposed 
into  two  components.  The  most  prevalent  is  a  diffuse  com¬ 
ponent.  This  has  instantaneous  intensity  that  typically  under¬ 
goes  random  fluctuations  that  obey  the  central  limit  theorem 
about  an  expected  value  that  decays  uniformly  with  time.  For 


^r(R,.,R0|r)=  Q(.f)®T(Rr,R0)e-i2”f'df,  (33) 

J  —00 

or  equivalently 

1R7-(R,,R0|f)=  J  J  J ''Pi(R).  — R|R0  — R|f)  dX  dY  dZ, 

V 

(34) 

where  ,'I/’S(R,.— R|R()  — R|f)  can  be  obtained  directly  from 
Eq.  (25)  for  distant  scatterers.  The  simplicity  of  this  equation 
is  deceptive.  While  it  can  be  evaluated  in  a  relatively 
straightforward  manner  for  deterministic  targets,  its  interpre¬ 
tation  and  implementation  become  far  more  difficult  for  sto¬ 
chastic  targets.  The  covariance 

<^r(R,,R0|f)^(R;,Rok')> 

-(^r(Rr,Ro|f))(^(R;,Ro|f')), 

for  example,  provides  a  second  moment  characterization  of 
the  reverberant  field  that  is  sufficient  for  most  remote  sensing 
applications,  where  a  statistical  correlation  between  scatter¬ 
ers  over  the  volume  V  is  implicitly  required  to  evaluate  the 
covariance.  This  quantity  is  most  useful  for  investigating  the 
performance  of  systems  that  employ  pulse  compression  in  an 
attempt  to  attain  high  temporal  resolution.  If  all  scatterers  are 
independent,  the  covariance  of  the  total  reverberant  field  in 
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time  can  be  obtained  directly  by  summing  the  covariances  of 
the  individual  scatterers  using  Eq.  (28). 

A  difficulty  arises,  however,  in  any  analytic  representa¬ 
tion  for  (^r(Rr,Ro|f)¥;|:(R;,Ro|f'))  when  expectation 
values  are  brought  within  the  modal  sums  and  spatial  and 
Fourier  integrals.  Even  for  a  single  random  scattering  patch, 
one  finds  that  the  second  moment  of  the  scatter-function, 
{S(0,4>\Si  ,4>i)\fS*{6'  ,4>' ',9'i  ,  0/ )  |/'),  is  required  for  joint 
correlation  across  both  wave  number  and  frequency.  We  are 
only  aware  of  analytic  derivations  existing  for  this  second 
moment  when  0=0' ,  </>=<// ,  6t=  0[ ,  r/>,  =  4>[ ,  and/=/',  as 
noted  in  the  definition  of  the  expected  bistatic  scattering 
cross  section  of  a  random  target  or  rough  surface  patch  in 
Appendix  A.  This  difficulty  is  circumvented  when 
(1T'r(R,.,R0|t)'vI,*(R' ,R0|t'))  is  estimated  by  Monte  Carlo 
simulation  since  'I''j(Rr,R0|f)''P*(R'  ,R0|f')  is  averaged 
across  independent  realizations  of  the  rough  surface.  Broad¬ 
band  scattering  statistics  can  then  be  obtained  by  either  the 
saddle  point  method  of  Sec.  IIE  or  the  Fourier  synthesis  of 
Eq.  (20)  since  only  the  deterministic  scatter  function  need  be 
known  to  compute  the  scattered  field  for  a  given  realization 
of  the  rough  surface.  In  this  way,  by  Monte  Carlo  simula¬ 
tions,  ('vPr(Rr,R0|f)''P*(R'  ,R0|r'))  can  be  readily  obtained 
using  the  theory  developed  in  Sec.  H  to  investigate  the  per¬ 
formance  of  systems  that  employ  high-resolution  temporal 
processing  such  as  pulse  compression. 

B.  Charting  diffuse  reverberation  when  system 
integration  time  spans  dominant  signal  energy 

A  simpler  analytic  approach  than  the  saddle  point  ap¬ 
proximation  or  Fourier  synthesis  can  be  employed  to  inves¬ 
tigate  system  performance  when  the  integration  time  of  the 
measurement  system  T  is  sufficiently  long  to  include  the 
dominant  signal  energy  returned  from  the  target  or  scattering 
patch.  In  this  case,  Parseval’s  theorem  can  be  applied  to  the 
Fourier  integral  of  Eq.  (20),  converted  to  absolute  coordi¬ 
nates,  to  obtain  the  time-averaged  mutual  intensity  expected 
at  R(.  from  a  target  or  scattering  patch  at  R  due  to  a  source  at 
Ko. 

1  ft+T/2 

/(R,R,.,R0,f)=  —  <|¥,(Rr-R|R0-R|r0)|2>*o 

1  Jt-T/2 

~  f  -R|R0-Rko)l2)  dt0 

=\ |°°Je(/)l2(|3>,(Rr-R|Ro-R)l2)^ 

(35) 

where  t—  772  is  less  than  or  equal  to  the  arrival  time  of  the 
scattered  signal. 

This  type  of  incoherent  integration  is  typically  used  in 
the  reception  of  narrow-band  source  waveforms,  and  is  also 
often  used  in  the  analysis  of  broadband  returns  from  explo¬ 
sive  sources  such  as  SUS22  where  the  exact  time  function  of 
the  source  is  unknown.  While  it  is  equally  valid  for  wave¬ 
forms  of  arbitrary  bandwidth,  it  does  not  take  advantage  of 


the  full  pulse  compression  possible  for  broadband  wave¬ 
forms.  However,  it  is  often  unclear  in  practice  whether  or  not 
pulse  compression  can  be  meaningfully  exploited  with  sig¬ 
nals  received  after  dispersive  waveguide  propagation.  For 
the  type  of  incoherent  time  average  specified  in  Eq.  (35),  the 
problems  mentioned  in  Sec.  Ill  A  are  alleviated  since  only 
(S(  6,  (f>\  0i ,  (/)/) | fS*(0' ,  4>' ;  0'i  ,  4>[)\ f)  need  be  evaluated  for 
a  given  random  scatterer,  so  that  Eqs.  (9)  and  (10)  can  be 
directly  applied,  since  frequency  cross  terms  vanish  as  a  re¬ 
sult  of  Parseval’s  theorem.  A  center  frequency  approxima¬ 
tion  to  Eq.  (35)  can  often  be  made  for  narrow-band  signals  as 
discussed  in  Appendix  B. 

If  the  time  spread  of  the  signal  due  to  dispersion  in  the 
waveguide  A  ts  is  small  compared  to  the  time  duration  Ts  of 
the  source  signal,  the  expected  horizontal  range  resolution 
Ap  of  the  system  will  take  roughly  the  same  form  as  in  free 
space  Ap  =  cTJ2,  for  narrow-band  signals,  where  c  is  the 
mean  horizontal  propagation  speed  of  the  signal  between 
source  and  receiver  in  the  waveguide.  In  this  case,  the  inte¬ 
gration  time  T  of  the  system  can  be  set  to  its  minimum  value 
of  Ts .  Both  A  ts  and  c  can  be  quantitatively  defined  in  terms 
of  the  received  field  as  in  Ref.  25.  They  depend  on  the  acous¬ 
tic  properties  of  the  waveguide,  the  signal  time  dependence, 
and  source-receiver  geometry.  For  the  narrow-band  ex¬ 
amples  of  Sec.  IV,  simulations  show  that  1500 m/s,  tc 
«»p  +  p0 ,  and  A  ts/Ts  is  small,  where  p  and  p0  are  defined  in 
Eq.  (30a). 

A  typical  bistatic  sonar  system  will  resolve  a  patch  of 
seafloor  A (R,R,. ,R0),  the  dimensions  of  which  depend  on 
the  receiving  array  aperture,  frequency,  and  the  bistatic  ge¬ 
ometry  of  the  source,  receiver,  and  seafloor  patch  as  dis¬ 
cussed  in  Appendix  C  of  Ref.  26  and  Refs.  27  and  28.  For  a 
monostatic  measurement  A  =  pApAcp,  where  A  cp=\/LA  is 
the  Rayleigh  resolution  of  the  horizontal  aperture  of  length 

La- 

For  convenience,  assume  that  the  horizontal  origin  in  an 
absolute  reference  frame  is  chosen  to  be  at  the  center  of  the 
receiving  array  Z,.=  (0,0 ,Zr).  Let  the  beamformed  output  of 
a  receiving  array  located  along  the  T,.-axis,  obtained  by  spa¬ 
tial  Fourier  transform  of  the  time-harmonic  scattered  field 
across  the  array  aperture,  be  denoted  by 


0>fi(<£s,Zr,R,R0)  =  f  r(7r)e(/)<&s(Rr-R|R0-R) 

J  —00 

Xeik  sin  <psrrdYr,  (36) 


where  i ps  is  the  azimuth  the  array  is  steered  towards,  <p  is  the 
azimuth  of  the  scattering  patch,  and  T(Yr)  is  the  array  taper 
function.  Suppose  a  uniform  rectangular  taper  function  is 
used  with  T(Yr)=  1  /LA  for  —LA/2^Yr^LA/2  and  zero  else¬ 
where,  and  the  seafloor  scattering  patch  is  in  the  far  field  of 
the  array,  such  that  |R|>LA/X,  and  the  scattering  patch  be¬ 
haves  as  a  point  target  to  the  array  so  that  the  angle  it  sub¬ 
tends  at  the  array  is  less  than  ~k/LA .  Under  these  assump¬ 
tions,  the  spectral  density  or  field  variance  received  from  this 
patch  can  be  well  approximated  by 
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Atmosphere 


cb%  db 3  ab3  Sediment  Half-Space 

FIG.  1.  The  geometry  of  the  waveguide  which  has  a  water  column  com¬ 
prised  of  upper  layer  sound  speed  cwX  for  0<Z<25,  lower  layer  sound 
speed  cw2  for  35<Z<100,  and  transition  layer  sound  speed  cwl  — (cwl 
—  cw2)(Z—  25)/10  for  25=SZ=S35.  The  water  column  density  is  dw 
=  1000  kg/m3  and  the  attenuation  is  arv„  =  6.0X  10-5  dB/X.  The  bottom  can 
have  up  to  two  sediment  layers.  The  upper  and  middle  sediment  layers  have 
respective  thicknesses,  sound  speeds,  densities,  and  attenuations  of  hl ,  cbl , 
dbl ,  abl  and  h2,  cb2,  db2,  ab2 ,  overlying  a  sediment  half-space  of  sound 
speed  cb 3 ,  density  db 3 ,  and  atteunation  a b2l .  The  monopole  source  is  colo¬ 
cated  with  receiving  array  center,  with  array  axis  normal  to  the  range-depth 
plane  of  the  sketch.  Source  and  receiver  may  be  placed  anywhere  in  the 
water  column.  The  submerged  target  may  be  placed  in  the  upper  or  lower 
layers  of  the  water  column  where  sound  speed  is  constant  as  indicated  in 
Fig.  2.  Seafloor  and  buried  riverbank  features  may  also  be  included  at  the 
water-sediment  and  sediment-layer  to  sediment  half-space  interfaces  as  in¬ 
dicated  in  Fig.  6.  Squiggly  lines  indicate  statistically  rough  interfaces. 


(|cDB(lp,Zf.,R,R0)|2)-|(<i>B(<p,Zr,R,R()))|2 

/  A  \2  oo  oo 

=  -H  |G(/)I2E  2  |Bm(Zr-R)|2 

\  k  J  m  =  1  n  =  1 

X  |A„(R0-R)|-Cmn(Zr-R,Zr-R|R0-R) 

sin{(LA/2)sin  <p[fc-Re{gm}]}  2 
(LA/2)sin<p|>-Re{£m}] 

upon  substitution  of  Eq.  (17)  into  Eq.  (36)  with  <f>s  =  cf>  so  the 
array  is  steered  toward  the  patch.  With  the  assumption  that 
the  resolution  footprint  A  is  much  larger  than  AA ,  the  area  of 
a  given  patch,  and  all  patches  are  statistically  independent, 
the  total  variance  of  the  received  field  from  seafloor  within 
the  system  resolution  footprint  can  be  written  as  the  sum  of 
the  variances  of  each  patch  via 


(a) 


Sediment  Half-Space 


(b) 


Sediment  Half-Space 


FIG.  2.  Three  scenarios  for  the  active  detection  of  a  submerged  pressure 
release  sphere  of  radius  a=  10  m.  The  water  column  is  modeled  as  either 
having  constant  sound  speed  or  as  downward  refracting.  The  bottom  is 
composed  of  either  a  pure  sediment  half-space  or  a  single  sediment  layer 
over  a  sediment  half-space,  (a)  Monopole  source  and  horizontal  receiving 
array  center  are  colocated  at  50-m  depth  with  target  at  50-m  depth  also,  (b) 
Source  and  receiving  array  center  are  colocated  at  10-m  depth  with  target  at 
50-m  depth,  (c)  Source  and  receiving  array  center  are  colocated  at  10-m 
depth  with  target  at  15-m  depth. 


When  Parse  val’s  theorem  is  invoked  again  under  the 
assumption  that  the  integration  time  of  the  measurement  sys¬ 
tem  includes  the  dominant  energy  returned  from  the  resolved 
patch  of  seafloor,  after  time-domain  beamforming  and  finite 
time  averaging  by  the  receiver  over  period  T,  the  field  vari¬ 
ance  from  seafloor  within  the  system  resolution  footprint  of 
area  A(R,R/.  ,R0)  centered  at  (X,Y)  becomes 

VB(X,y)=^f  VB{XJ)df.  (39) 


Vb[X,Y)=  f  f  «|<f>B(<p,Z,.,R',R0)|2) 

A(R,Zr,Ro) 

-\<‘PB(cp,Zr,R',R0))\2)^XjdX'  dY', 

(38) 

where  n  =  (nx,ny  ,nz)  is  the  surface  normal  at  R.  Since  the 
differential  area  dX' dY'  must  be  normalized  by  the  horizon¬ 
tal  projected  area  of  each  potentially  inclined  patch  at  R'  to 
allow  horizontal  integration,  Eq.  (38)  does  not  allow  vertical 
patches. 


Equations  (37)-(39)  imply  a  reduction  in  reverberation  level 
for  off-broadside  beams  due  solely  to  modal  dispersion.  Only 
at  broadside  does  the  phase  speed  of  the  incident  waves 
match  that  expected  under  the  nondispersive  assumptions  of 
plane  wave  beamforming.  Only  broadside  beamforming  is 
considered  in  the  simulations  of  the  present  article  to  elimi¬ 
nate  this  effect  from  the  analysis.  The  effects  of  modal  dis¬ 
persion  on  beamformed  reverberation  are  investigated  in 
Ref.  25. 

Under  the  present  assumptions,  reverberation  measured 
in  time  can  be  charted  in  space  for  any  bistatic  geometry 
using  a  look-up  table  comprised  of  the  mean  time  delay  from 
source  to  scattering  patch  r(r,  ,r0)  and  scattering  patch  to 
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receiver  r(r,rf).25  Similarly,  reverberation  modeled  with  the 
spatial  formulation  of  Eqs.  (38)  and  (39)  can  be  made  a 
function  of  time  by  reversing  the  procedure. 

C.  Lambert-Mackenzie  scattering  and  reverberation 

A  number  of  simplifications  are  possible  when  the  scat¬ 
tering  surface  has  Lambertian 29  behavior.  The  covariance  of 
the  scattering  function  for  a  Lambertian  scattering  patch  of 
area  AA  with  albedo  e  takes  the  form  of  Eq.  (12a)  with 

c  Zm(Rr  -  R,  R;  -  R|  R0  -  R)  =  k2- 1  isl ■  n\\  iim  -n\AA, 

(40) 

where  iim  and  ;  v/  are  the  directions  in  which  the  down-going 
component  of  the  mth  incident  and  up-going  component  of 
the  /th  scattered  modes  propagate  and  h  is  the  seafloor  nor¬ 
mal,  pointing  away  from  the  water  column,  in  a  scattering- 
patch-centered  coordinate  system  where  the  positive-z  axis 
points  downward.  The  differential  scattering  cross  section  of 
the  surface  patch  is  then  given  by  the  product  of  4  it  Ik 2  and 
the  right-hand  side  of  Eq.  (40). 

Under  Lambertian  scattering,  the  cross  spectral  density 
of  the  scattered  field  given  in  Eq.  (17)  then  becomes  express¬ 
ible  in  terms  of  single  summations 

<3>,(r)3>*(r'))-<®,(r))<®*(r')} 

S,(r)B?(r')|i,rn| 


=  167reAA 


2 

/=] 


The  surface  projection  factors  can  be  written  in  terms  of  the 
incident  and  scattered  wave  number  components,  the  mea¬ 
surement  geometry,  and  the  orientation  of  the  surface  patch 
via 


\isrn\ 


6r  .  ■  I,7' 

—  [n^cos  <p  +  ny  sin  (p J+  —  nz 


(42a) 


%  [nx  cos(  cp0  +  tt)  +  ny  sin(  <p0  +  tt)  ]  -  -2  nz 


(42b) 


The  scattered  field  covariance  from  a  given  seafloor  patch  in 
a  target-oriented  frame  is 


(<&i(r|r0)<hf  (r'|r0))  — (<J>J(r|r0))(<i)f  (r'|r0)) 


=  167reAA  2  B,(r)Bf(r') 


1=  l 


—  [nx  cos  (p  +  ny  sin  <p] 


,  Vi 
+  J7'h 


2  |Am(r0)|2 

m=  1 


y  [«aCOs(<Po+7t) 


+  nysin(<p0+7r)]-  -ynz 


(43) 


X  2  l^m(r0)|2|tjm-n 

\m  =  1 


(41) 


where  nx=ny=  0  for  a  bottom  with  zero  mean  inclination.  In 
an  absolute  frame  it  becomes 


($s(Rr-R|R0-R)<Ds*(R;-R|R0-R))-($J(Rr-R|R0-R))($f(R;-R|R0-R)) 


=  167reAA 


2  R/(Rr-R)Bf(R;-R) 

1=  1 


g,nx(Xr-X)  +  ny{Y-Y)  |  y,  ^ 
k  j(X-X)2+(Y-Y)2  +  /cHz  ) 


X 


oc 


2 

m  =  1 


K,(r„-r)|2 


g,„  nx(XQ-X)  +  ny(Y0-Y)  _  ym  ^ 
k  ^(Xq-X)2  +  (Y0-Y)2  k  ) 


(44) 


where  the  components  of  the  surface  normal  n  are  now  a  function  of  the  X,  Y,  Z  position  of  the  surface  patch  center.  Under 
far-field  assumption,  the  field  variance  received  from  seafloor  within  the  system  resolution  footprint  centered  at  ( X,Y )  and 
averaged  over  time  period  T  can  be  well  approximated  by 

6  nx(X-X')  +  ny(Y-Y')  |  y, 
k  nz^l(Xr-X')2  +  (Yr-Y')2  k 

A(R,Zr  ,Rq) 


VB{X,Y)=l^j°°jQ(f)  |2  fj 


2  |B,(Z,.-R')|: 

1=1 


sin((LA/2)sin  ^[k-Relg,}])  2\ 

(La/2)  simp[k-Re{^,}]  J 


00 


2 

m  =  1 


|A,„(R0-R')|2 


gm  r,x(X0-X')  +  ny(Y0-Y')  ^  ym 
k  nzJ(X0-X')2  +  (Y0-Y')2  k 


dX'  dY'  df. 


(45) 


918  J.  Acoust.  Soc.  Am.,  Vol.  109,  No.  3,  March  2001  N.  C.  Makris  and  P.  Ratilal:  Unified  model  for  reverberation  and  scattering  918 


Sand,  Sphere  Backscatter 
Silt,  Sphere  Backscatter 
Sand,  Seafloor  Reverberation 
Silt,  Seafloor  Reverberation 


1 1lk  H  fl  /!  A  \  i  t  , 

mJkMWUAj 


10  15  20  25  30  35 

Range  (km) 


10  15  20  25  30  35  40  45  50 

Range  (km) 


°3-100 


.J 

T3 


ft. 

■o 


-  2m  Sand  over  Silt,  Sphere  Backscatter 

-  5m  Sand  over  Silt,  Sphere  Backscatter 

-  -  -  2m  Sand  over  Silt,  Seafloor  Reverberation 
—  5m  Sand  over  Silt,  Seafloor  Reverberation 


1.2  5  10  15  20  25  30  35  40  45  50 

Range  (km) 


FIG.  3.  The  scattered  field  from  a  submerged  pressure-release  sphere  of  radius  a  =  10  m,  at  /=  300  Hz  and  center  at  50-m  depth,  and  Lambert-Mackenzie 
reverberation  from  the  seafloor  within  the  broadside  resolution  footprint  of  the  monostatic  system  as  a  function  of  range  for  a  water  column  with  constant 
sound  speed  of  1500  m/s,  i.e.,  cwi=cw2  =  1500  m/s.  Monopole  source  and  receiving  array  center  are  colocated  at  50-m  depth.  Range  increases  along  the  .v-axis 
and  depth  along  the  z-axis,  with  the  array  axis  along  the  y-axis.  Source  strength  is  0  dB  re  1  /jPa  @  1  m.  Reverb  modeled  with  /  1/2  s  duration  cw  source 

signal  at  300  Hz  and  receiving  array  resolution  XI L  =  3.7  degrees,  (a)  Pekeris  waveguide  examples  for  bottom  half-spaces  composed  of  either  sand  or  silt,  i.e., 
h1=h2  =  0.  (b)  Bottom  has  a  silt  layer  of  either  li{  =  2  m  or  hl  =  5  m  overlying  a  sand  half-space,  and  h2  =  0.  (c)  Bottom  has  a  sand  layer  of  either  h j 
=  2mor/i1  =  5m  overlying  a  silt  half-space,  and  h2  =  0.  Error  bars  show  the  5.6  dB  standard  deviation  in  reverb  level. 


when  T  is  sufficiently  large  for  the  dominant  energy  of  the 
scattered  field  to  be  received.  This  result  for  the  Lambertian 
seafloor  offers  significant  advantages  in  implementation 
through  the  separation  of  the  incident  and  scattered  modal 
summations.  For  narrow-band  waveforms,  terms  within  the 
modal  summations  of  Eq.  (45)  often  vary  so  slowly  that  they 
can  be  approximated  as  a  constant  function  of  frequency 
over  the  dominant  portion  of  the  spectral  window  Q{f).  This 
greatly  simplifies  computations  as  shown  in  Appendix  B. 


Let  the  x  and  y  components  of  the  gradient  of  the  surface 
Zj(x,y)  be  denoted  by 

dzs 

P=~,  (46a) 

q=j1-  (46b) 

dy 

The  surface  normal  can  be  expressed  as 


D.  Perturbation  theory  for  diffuse  rough  surface 
scattering  and  reverberation 

Perturbation  theory  can  also  be  used  to  calculate  the 
field  scattered  by  a  rough  surface.  The  advantage  of  pertur¬ 
bation  theory,  when  it  is  applicable,  is  that  it  is  derived  from 
first  principles  and  so  requires  knowledge  of  only  the  geo¬ 
acoustic  properties  of  the  media,  such  as  sound  speed  and 
density,  as  well  as  a  second  moment  characterization  of  the 
statistical  properties  of  the  scattering  surface. 


_  (-p,-qA) 

V 1  +p2  +  q 2 

along  with  two  orthonormal  surface  tangents 

_  (i,o,p) 
fl  VTV’ 

and 


(47a) 


(47b) 
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FIG.  4.  Same  as  Fig.  3  except  water  column  is  layered  with  c w |  =  1520  m/s  and  cw2=  1500  m/s  and  monostatic  source-receiver  as  well  as  target  sphere  are 
at  variable  depth.  Only  the  cases  of  pure  sand  or  pure  silt  bottom  half-spaces  are  shown,  (a)  Source-receiver  and  sphere  center  are  at  50-m  depth,  (b) 
Source-receiver  are  at  10-m  depth  while  sphere  center  is  at  50-m  depth,  (c)  Source-receiver  are  at  10-m  depth  while  sphere  center  is  at  15-m  depth. 


(- pq,\+p2,q ) 
t2~^(\+p2)(l+p2  +  q2)’ 


(47c) 


where  is  tl  obtained  by  taking  an  infinitesimal  step  along  the 
surface  on  the  x  axis,  and  t2  is  the  crossproduct  of  n  and  tx . 

The  projections  of  the  incident  and  scattered  wave  num¬ 
ber  vectors  on  the  surface  then  become 


\K\2=(k-tl)2  +  (k-t2)2.  (51b) 

A  plane  wave  incident  from  medium  1  half-space  that  is 
reflected  from  strata  below  has  total  reflection  coefficient30 


r(K,.)= 


r12(K/)+r'(K,-)ef2y(2>*i 

l  +  rI2(Ki)r'(K1.)e'‘2^2,*i’ 


(52a) 


Ki  =  (krtl)tl  +  (krt2)t2, 

(48a) 

where 

K=(k-t1)tl  +  (k-t2)t2, 

(48b) 

p2l^k\-K2i-pxl^k1x-K2i 

where,  for  incident  mode  n  and  scattered  mode  m. 

12(K'j  p2l4k\-K2  +  pxl4k\-K2 

(52b) 


*;  =  (£„  cos{cf>0+  ir),£„  sin(</>0+  tt),  yn),  (49a) 

k  =  ( cos  </>,  sin  <f>,  ym ) ,  (49b) 

so  that 

K;.  •  K=  (krtMk- tl)  +  (krt2)(k-t2),  (50) 

and,  for  example, 

|K1-|2  =  (*i-*1)2  +  (*,-»2)2.  (51a) 


is  the  reflection  coefficient  from  the  medium  1  to  medium  2 
interface,  r'(K,)  is  the  total  reflection  coefficient  from  all 
strata  below  medium  2  for  a  plane  wave  incident  from  me¬ 
dium  2,  li  [  is  the  thickness  of  the  layer  containing  medium 
2, and  y(21  is  the  vertical  wave  number  component  of  me¬ 
dium  2. 

The  differential  scattering  cross  section  of  a  surface 
patch  of  area  AA,  from  first-order  perturbation  theory,  can 
be  expressed  as30 
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FIG.  5.  Pekeris  waveguide  with  varying  source-receiver  and  target  depth,  (a)  Same  as  Fig.  3(a)  except  source-receiver  is  at  10-m  depth  and  target  is  at  50-rn 
depth,  (b)  Same  as  Fig.  3(a)  except  source-receiver  is  at  10-m  depth  and  target  is  at  15-m  depth. 


fFpertC  ot,  /?;  «/,/?,■)  =  4  IT  A  A 


|[r(K)+i][r(K,)+i]|2 


K-K;  P]{  k,k,-,n)\ 
k\  dbt  I 


2 


W(K-Kj), 


(53) 


where 

r2(KK)-d2  |k;'n|  lk  nl  ( 1  —  ^(K')\ [  i  ~r(K)\ 
bt  kt  k  \i+r(Ki)/\i+r(K)/’ 

(54) 

and 

db,=db/dt,  (55) 

K—  Kt  I  Kb  ,  (56) 

where  d,  and  d h  are  the  respective  densities  above  and  below 
the  scattering  interface  and  k,  and  Kh  the  respective  wave 
number  magnitudes  above  and  below  the  scattering  interface. 

Following  Moe  and  Jackson,30  the  roughness  of  the 
given  surface  patch  is  assumed  to  follow  the  isotropic  power 
law 

JF(K)  =  w2|K|~y.  (57) 

With  the  assumption  that  the  scattering  patch  is  much  greater 
than  the  wavelength  so  that  the  incident  and  scattered  modes 
are  decorrelated  by  the  scattering  process,  the  covariance  of 
the  scatter  function  is  given  by  Eq.  (12a)  with 


C  mn(  U  l®*o)  a  tTpert(  Ci.m  a  &  *  ,  (f>Q  T  tt)  .  (58) 

47T 

Upon  substituting  Eq.  (58)  into  Eqs.  (17),  or  (38)  after  beam¬ 


forming,  it  is  found  that  the  covariance  of  the  field  scattered- 
from  a  rough  surface  patch  that  obeys  first-order  perturbation 
theory  involves  a  double  summation  over  the  waveguide 
modes.  Evaluating  this  is  significantly  more  computationally 
intensive  than  the  product  of  single  modal  summations  found 
in  the  Lambert-Mackenzie  formulation. 

E.  Coherent  reverberation  from  deterministic  and 
stochastic  geological  features 

There  are  two  general  kinds  of  seafloor  scatterers  that  do 
not  decorrelate  the  incident  or  scattered  modes.  A  seafloor 
scatterer  of  the  first  kind  can  be  modeled  as  a  deterministic 
feature,  with  known  or  computable  far-field  scatter  function, 
that  can  have  arbitrary  size  compared  to  the  wavelength  so 
long  as  it  falls  within  the  resolution  footprint  of  the  active 
sonar  system.  The  feature  must  be  distinct  from  the  other¬ 
wise  range-independent  boundaries  of  the  stratified  medium 
in  order  to  induce  scattering. 

A  compelling  canonical  example  of  a  seafloor  scatterer 
of  the  first  kind  is  a  smooth  flat  inclined  segment  of  the 
seafloor,  such  as  a  seafloor  or  subseafloor  river  channel,  ice¬ 
berg  scour,  or  submerged  hillside,  that  can  be  modeled  as  a 
flat  plate  with  scattering  characteristics  determined  by  its 
size,  inclination,  and  the  local  geo-acoustic  properties  of  the 
interface.  The  3-D  scatter  function  for  a  rectangular  surface 
patch  with  total  reflection  coefficient  r(K,),  for  example, 
can  be  readily  determined  by  applying  Green’s  theorem,  Eq. 
(Al),  for  a  plane  wave,  with  wave  number  magnitude  k\ , 
incident  in  the  direction  ( «,■  ,/?,)  and  a  far-field  receiver  in 
the  direction  (a,  0)  with  respect  to  the  patch  centroid.  If  the 
patch  is  assumed  to  be  at  inclination  \  from  horizontal, 
where  the  angle  x  comprises  a  counter-clockwise  rotation 
about  the  y  axis,  the  scatter  function  takes  the  form 
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FIG.  6.  Scenarios  for  the  active  detection  of  seafloor  and  subseafloor  riverbank  features.  The  water  column  is  modeled  as  having  constant  sound  speed,  i.e., 
cwi  =  cw2=  1500  m/s,  and  monopole  source  and  horizontal  receiving  array  center  are  colocated  at  50-m  depth  in  all  cases,  (a)  Bottom  is  sediment  half-space 
with  seafloor  feature,  i.e.,  hl=h2= 0.  (b)  Bottom  is  composed  of  a  single  sediment  layer  with  double-interface  seafloor  feature,  i.e.,  h2  =  0.  (c)  Bottom  is 
composed  of  a  single  sediment  layer  with  subseafloor  feature,  i.e.,  h2  =  0.  (d)  Bottom  is  composed  of  two  sediment  layers  with  double-interface  subseafloor 
feature. 


S(a,p-,ai,0i)  = 


.  *1 
i  - — 
47 T 


LxLy [( 1  —  r(K,)){cos  ctj  cos  x+  sin  or,  sin  x  cos  /?,■}  +  (!  +T(K()){cos  a  cos  sin  a  sin  x  cos  /3}] 


k \LX  r  .  . 

^1  Ly  . 

X  sine 

— - —  { ( sin  a(  cos  /3;  cos  x  ~  cos  or,  sin  —  ( sin  a  cos  / 3  cos  x  ~  cos  a  sin  x) } 

sine 

^  {sin  or,-  sin  f3j  —  sin  a  sin  f3\ 

(59) 


where  sinc(x)  is  defined  as  sinx/x.  The  reflection  coefficient 
can  be  determined  from  Eq.  (52)  with  the  understanding  that, 
in  the  present  geometry,  the  squared  magnitude  of  the  trans¬ 
verse  component  of  the  incident  wave  number  vector  on  the 
inclined  surface  patch  is 

Kj  =  k\  —  k2(cos  a,  cos  sin  a,  sin  \  cos  /3,- ) 2 .  (60) 

Irregularities  in  the  surface  can  make  its  scatter  function  de¬ 
viate  from  that  given  in  Eq.  (59).  For  realistic  seafloor  and 
subseafloor  riverbanks,  however,  it  is  reasonable  to  assume 
that,  for  a  low-frequency  active  system28  at  long  range  in  a 
shallow  water  waveguide  where  propagation  is  near  horizon¬ 
tal,  the  product  of  the  amplitude  of  such  irregularities  and  the 
normal  component  of  the  wave  number  vector  with  respect 
to  the  surface  will  be  small  enough  that  the  irregularities  will 
have  a  negligible  effect  on  the  field  from  the  riverbank. 


A  seafloor  scatterer  of  the  second  kind  is  a  randomly 
rough  rather  than  deterministic  feature  but  is  appropriately 
modeled  with  completely  coherent  modes  when  the  ratio  of 
wavelength  to  system  range  resolution,  X/Ap,  is  near  or 
greater  than  the  equivalent  vertical  propagation  angle  of  the 
highest  order  trapped  mode,  which  in  many  shallow  water 
scenarios  is  roughly  the  bottom  critical  angle.  This  situation 
occurs  for  active  sonar  systems  with  high  range-resolution 
and  can  lead  to  the  formation  of  range-dependent  rings  in 
charted  reverberant  intensity  caused  by  modal  interference.24 
(Lepage7  has  recently  described  scenarios  in  which  such 
rings  can  form  even  in  narrow-band  reverberation  at  short 
ranges.)  The  level  of  returns  can  be  estimated  by  appropri¬ 
ately  modeling  the  seafloor  scatter  function.  If  the  system 
resolution  footprint  extends  over  many  wavelengths  in  any 
direction  and  the  correlation  length  of  surface  roughness  is 
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much  smaller  than  the  system  resolution  footprint,  then  the 
scattering  function  for  the  seafloor  over  this  area  can  be 
treated  as  a  fluctuating  target.  If  the  resolution  footprint  is  on 
the  order  of  the  wavelength  or  the  correlation  length  of  sur¬ 
face  roughness,  a  quasi-deterministic  description  of  the  scat¬ 
tering  process  can  be  used. 

IV.  ILLUSTRATIVE  EXAMPLES 

In  all  the  illustrative  examples  of  this  section,  a  water 
column  of  100-m  depth  is  used  to  simulate  a  typical  conti¬ 
nental  shelf  environment.  The  sound  speed  structure  of  the 
water  column  varies  from  iso-velocity  to  downward  refract¬ 
ing  layers  with  constant  density  of  1  g/cm3  and  attenuation  of 
6. OX  10~5  dB/A.  The  seabed  is  comprised  of  sand  or  silt 
half-spaces,  with  up  to  two  sediment  layers,  comprised  of 
sand  or  silt,  over  a  sand  or  silt  half-space.  The  density,  sound 
speed  and  attenuation  are  taken  to  be  1.9  g/cm3,  1700  m/s, 
and  0.8  dB/A  for  sand,  1.4  g/cm3,  1520  m/s,  and  0.3  dB/A  for 
silt,  and  1.2  g/cm3,  1510  m/s,  and  0.3  dB/A  for  light-silt. 
Scattering  and  reverberation  calculations  are  made  for  a  sub¬ 
merged  target,  roughness  at  the  water-seabed  interface, 
roughness  at  the  interface  between  the  upper  seabed  layer 
and  lower  half-space,  as  well  as  for  anomalous  features  of 
the  seafloor  or  subseafloor  that  return  geological  clutter.  The 
latter  are  taken  to  be  seafloor  river  banks  at  the  water-seabed 
interface  or  subseafloor  riverbanks  at  the  interface  between 
the  upper  sediment  layer  and  lower  half-space.  The  geometry 
of  the  waveguide  is  sketched  in  Fig.  1. 

A  horizontal  line  array  with  N=  32  equally  spaced  ele¬ 
ments  of  length  La  =  (N—\)\I2  at /=  300 Hz  is  used  as  a 
receiver  and  a  cw  pulse  of  T=\-  s  duration  centered  at  / 


FIG.  7.  The  field  at  /=  300  Hz  scattered  from  a  coherently  scattering  rect¬ 
angular  patch  of  area  100X100  m2  representing  a  seafloor  riverbank  for 
scenario  shown  in  Fig.  6(a),  constant  sound  speed  water  column  over  pure 
silt  or  sand  half-spaces.  Range  increases  along  the  x  axis  and  depth  along  the 
z  axis.  The  square  riverbank  surface  has  two  edges  parallel  to  the  y  axis,  and 
is  inclined  10°  from  the  x  axis.  Constant  sound  speed  in  the  water  column  is 
assumed  for  all  examples  with  cwl  =  cw2=  1500  m/s.  Lambert-Mackenzie 
reverberation  within  the  range-dependent  resolution  footprint  of  the  mono¬ 
static  system  is  also  shown  separately  for  the  water- sediment  interface  (sea¬ 
floor).  Source  strength  is  0  dB  re  1  /jlP a  @  1  m.  Diffuse  reverb  modeled 
with  T=  1/2  s  duration  cw  source  signal  at  300  Hz  and  receiving  array 
resolution  A./L  =  3.7  degrees. 


=  300  Hz  is  used  as  a  source  waveform  for  all  simulations  of 
diffuse  reverberation.  Targets  beyond  LAI A  are  in  the  far 
field  of  the  array,  which  begins  at  roughly  1.2  km.  The 
beamformed  field  from  an  object  that  falls  within  the  broad¬ 
side  beam  of  the  array,  in  the  absence  of  other  sources  or 
scatterers,  equals  the  field  received  from  that  object  by  a 
single  hydrophone  at  the  array  center  when  Eq.  (36)  is  used 
with  uniform  taper  T(Yr)  =  l/LA.  If  the  same  object  is 
placed  at  the  same  range  but  within  an  off-broadside  beam,  a 
reduction  in  the  beamformed  output  may  occur  due  to  modal 
dispersion,  as  is  discussed  in  detail  in  Ref.  25.  For  simplicity, 
only  objects  and  reverberation  within  the  broadside  beam  are 
considered  in  the  present  article.  Only  monostatic  scenarios 
are  considered,  where  the  source  is  located  at  the  center  of 
the  receiving  array.  This  leads  to  a  range-dependent  resolu¬ 
tion  footprint  A  =  pApA<p,  where  c/p  =  c772=375m  and 
dtp=  A/L^ 3.7  degrees  for  the  given  array,  frequency,  and 
cw  pulse  length. 

A  center  frequency  approximation,  at  /=300Hz,  is 
made  for  all  scattering  calculations.  For  reverberation  calcu¬ 
lations  this  approximation  differs  from  the  full  spectral  inte¬ 
gration  by  less  than  0. 1  dB  for  the  examples  shown.  As  may 
be  expected  in  coherent  scattering  from  targets  where  modal 
interference  is  significant,  some  range-dependent  nulls  and 
valleys  in  the  sound  pressure  level  of  the  received  field  found 
in  the  single  frequency  calculation  may  be  partially  filled 
when  the  full  bandwidth  is  used  for  the  narrow-band  wave¬ 
forms  considered.  Since  this  filling  is  window  dependent,  as 
shown  in  Appendix  B,  only  center  frequency  calculations  are 
presented  in  the  main  text.  It  is  also  shown  in  Appendix  B 
that  in  some  valleys  of  some  single  frequency  calculations 
the  target  returns  may  fall  below  the  expected  reverberation 
level  but  will  be  above  this  level  when  the  full  bandwidth  of 
a  given  narrow-band  window  function  is  employed. 

Only  the  empirical  Lambert-Mackenzie  model  is  used 
in  comparisons  between  seafloor  reverberation  and 
submerged-object  returns  since  insufficient  data  on  the  req¬ 
uisite  environmental  parameters  at  low  frequency  are  avail¬ 
able  to  make  a  similar  comparison  with  perturbation  theory 
meaningful.  Perturbation  theory  calculations  are  only  used 
self-consistently  to  make  inferences  about  the  relative  level 
of  returns  from  different  kinds  of  seafloor  scatterers. 

A.  Submerged  target  echo  versus  diffuse 
reverberation  level  for  varying  source-receiver 
depth,  target  depth,  water  column,  and 
bottom  stratification 

The  geometry  for  active  detection  of  a  sphere  sub¬ 
merged  in  an  ocean  waveguide  is  sketched  in  Fig.  2  for  the 
illustrative  examples  of  this  section.  The  geometry  is  mono¬ 
static  with  co-located  omni-directional  point  source  and  re¬ 
ceiving  array  centers  at  50-m  depth.  The  sphere  center  is  also 
at  D  =  50-m  depth  at  array  broadside  with  variable  horizontal 
range.  The  field  back  scattered  from  a  pressure  release  sphere 
of  radius  a=10m  at  /=  300  Hz  is  shown  as  a  function  of 
range  in  Figs.  3(a)-(c)  in  decibels,  i.e.,  201og|<I>J|,  for  vari¬ 
ous  bottom  types  under  a  water  column  with  constant  sound 
speed  cH,=  1500m/s,  where  cwl  =cw2=  1500m/s.  The  scat¬ 
tered  field  is  computed  by  Eq.  (1),  with  scatter  function 
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FIG.  8.  Same  as  Fig.  7  except  single  layer  bottom  scenarios  of  Figs.  6(b)  and  (c)  for  coherent  seafloor  and  subseafloor  riverbank  scattering  are  investigated. 
Lambert-Mackenzie  reverb  within  the  sonar  resolution  footprint  is  also  shown  for  the  sediment-layer  to  sediment  half-space  interface  (subseafloor),  (a) 
Seafloor  riverbank  with  the  upper  sediment  layer  composed  of  silt  with  h{  —  2  m  and  the  lower  sediment  half-space  composed  of  sand.  Coherent  riverbank 
scattering  is  from  the  double  interface  of  water  to  silt  to  sand,  (b)  Subseafloor  riverbank  with  the  upper  sediment  layer  composed  of  silt  with  hx  =  2  m  and  the 
lower  sediment  half-space  composed  of  sand.  Coherent  riverbank  scattering  is  from  the  single  silt  to  sand  interface,  (c)  Seafloor  and  subseafloor  riverbanks 
as  in  (a)  and  (b)  but  with  the  upper  sediment  layer  now  at  hx  =  5  m  thickness,  (d)  Same  as  (c)  but  with  the  upper  sediment  layer  composed  of  sand  with 
h !  =  2  m  and  the  lower  sediment  half-space  composed  of  silt,  (e)  Same  as  (d)  except  h{  =  5  m. 


given  by  Eqs.  (8)  and  (9)  of  Ref.  11  with /(n)  replaced  by 
( —  1  )"/(n)  to  convert  from  Ingenito’s  definition  to  the  stan¬ 
dard  one  described  in  Sec.  II  A.31 

The  variance  of  the  field  scattered  from  the  seafloor 


within  the  range-dependent  resolution  footprint  of  the  sonar 
system  under  the  Lambert-Mackenzie  assumption  of  Eq. 
(45)  is  also  shown  in  Fig.  3  in  decibels,  i.e.,  101ogVB. 
Modal  interference  is  absent  due  to  the  modal  decoupling 
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assumed  in  diffuse  scattering  from  large  seafloor  patches. 
Ambiguous  returns  from  both  sides  of  the  line  array  are  in¬ 
cluded. 

The  scattered  field  from  both  the  target  sphere  and  sea¬ 
floor  is  highly  dependent  on  the  geo-acoustic  parameters  of 
the  bottom,  as  is  evident  in  Fig.  3(a)  where  significant  dif¬ 
ferences  arise  when  the  bottom  type  is  changed  from  sand  to 
silt.  The  differences  arise  primarily  because  the  number  of 
trapped  modes  is  significantly  larger  for  the  sand  half-space 
due  to  the  higher  critical  angle  of  28.1°  for  water  to  sand  as 
compared  with  the  9.3°  for  water  to  silt.  This  leads  to  a 
correspondingly  higher  mean  level,  of  roughly  20  dB,  for 
both  target  and  seafloor  backscatter  and  a  shorter  modal  in¬ 
terference  length  scale  in  the  scattered  field  from  the  sphere. 

In  the  Pekeris  waveguide  examples  of  Fig.  3(a),  the  tar¬ 
get  stands  tens  of  decibels  above  the  expected  reverberation 
within  the  broadside  beam  regardless  of  whether  the  bottom 
is  composed  of  sand  or  silt.  This  signal  excess  is  well  above 
the  reverberation  level  standard  deviation  of  5.6  dB  assum¬ 
ing  the  seafloor  scattering  obeys  circular  complex  Gaussian 
statistics,  in  accord  with  the  central  limit  theorem.21  If  a 
single  omni-directional  hydrophone  placed  at  the  center  of 
the  receiving  array  replaces  the  full  array,  the  reverberation 
levels  are  augmented  by  roughly  1 0 1  og(2  Tfld cp)  ~~  20  d  B  in 
Fig.  3.  The  target  sphere  then  no  longer  consistently  stands 
above  the  expected  reverberation  even  at  short  ranges,  for 
example,  within  a  few  kilometers.  A  directional  array  is  then 
necessary  to  spatially  filter  the  target  from  omni-directional 
reverberation  so  that  detection  can  be  practically  achieved  in 
the  given  scenarios. 

The  effect  of  bottom  properties  on  both  submerged  tar¬ 
get  scattering  and  reverberation  is  again  evident  when  lay¬ 
ered  bottoms  are  considered.  For  the  silt-over-sand  scenarios 
of  Fig.  3(b),  the  characteristics  of  the  field  scattered  from  the 
target  are  a  combination  of  those  found  for  the  silt  and  sand 
half-spaces.  As  the  silt  layer  increases  from  roughly  one-half 
to  a  full  wavelength,  the  rate  of  modal  interference  de¬ 
creases,  as  does  the  overall  level  of  the  scattered  field  from 
both  the  target  and  bottom.  When  the  layer  thickness  reaches 
a  full  wavelength,  the  level  of  reverberation  approaches  that 
obtained  for  a  pure  silt  bottom  as  range  increases.  The  low 
critical  angle  between  the  water-silt  interface  enables  greater 
bottom  penetration  than  is  possible  with  a  water-sand  inter¬ 
face.  The  high  attenuation  of  the  silt  layer  then  leads  to  bot¬ 
tom  loss  that  increases  with  the  thickness  of  the  layer. 

For  the  sand-over-silt  scenarios  of  Fig.  3(c),  the  field 
scattered  from  the  target  greatly  resembles  that  obtained  for 
the  pure  sand  bottom  of  Fig.  3(a).  The  match  becomes  better 
as  the  sand  layer  increases  in  thickness  from  one-half  to  a 
full  wavelength,  in  which  case  the  reverberation  increases 
from  a  few  decibels  below  to  roughly  the  level  found  for  a 
pure  sand  bottom.  In  the  latter  case,  the  silt  half-space  is 
effectively  insulated  from  the  water  column  by  evanescent 
decay  of  the  trapped  modes  in  the  sand  layer. 

The  absolute  and  relative  levels  of  target  and  reverbera¬ 
tion  echo  returns  are  highly  dependent  upon  the  water  col¬ 
umn  sound  speed  structure  as  well  as  source,  receiver,  and 
target  depth.  To  illustrate  this,  consider  the  typical  shallow 
water  downward  refracting  profile  shown  in  Fig.  1,  with 


cwi=  1520m/s,  cW2=1500m/s,  and  a  linear  transition  re¬ 
gion  in  between,  that  is  similar  to  what  is  found  in  continen¬ 
tal  shelf  waters  in  late  spring  and  summer  months.  Mono¬ 
static  measurements  of  the  field  scattered  from  a  10-m-radius 
pressure-release  sphere  are  again  made  with  the  same  array 
and  cw  tone  used  in  the  previous  examples.  The  target  is  at 
array  broadside  and  both  target  returns  and  reverb  within  the 
broadside  beam  are  plotted  as  a  function  of  range,  where  the 
reverb  is  computed  again  for  a  T=\ s  cw  at  300  Hz  center 
frequency.  Three  combinations  of  monostatic  source- 
receiver  and  target  depths  are  considered,  as  illustrated  in 
Fig.  2. 

First  consider  the  case  in  Fig.  4(a),  where  the  source, 
receiving  array,  and  sphere  center  are  at  50-m  depth  just  as  in 
the  Pekeris  waveguide  examples.  For  the  sand  bottom,  the 
levels  are  similar  to  those  found  in  the  corresponding  ex¬ 
ample  of  Fig.  3(a),  with  the  target  standing  out  by  tens  of 
decibels.  For  the  silt  bottom,  the  target  still  stands  tens  of 
decibels  above  the  reverberation  but  the  absolute  levels  of 
the  scattered  fields  decay  more  rapidly  with  range  in  the 
present  scenario  since  the  downward-refracting  profile 
causes  more  acoustic  energy  to  penetrate  into  the  bottom. 

Loss  of  energy  to  the  bottom  is  augmented  when  the 
source  and  receiver  array  are  placed  in  the  mixed  layer,  at 
10-m  depth,  while  the  target  remains  with  center  at  50-m 
depth,  as  shown  in  Fig.  4(b).  The  absolute  levels  of  both  the 
field  scattered  from  the  target  sphere  and  the  seafloor  are 
reduced  by  tens  of  decibels  beyond  a  few  kilometers’  range 
for  the  silt  bottom.  For  the  sand  bottom,  the  reverberation 
level  is  not  significantly  changed  by  moving  the  source  and 
receiver  into  the  mixed  layer.  Returns  from  the  target  sphere, 
however,  no  longer  stand  prominently  enough  above  the  ex¬ 
pected  reverberation  to  insure  detection,  given  a  5.6-dB  stan¬ 
dard  deviation  in  reverberation  level. 

The  situation  for  detection  again  changes  when  the  tar¬ 
get  sphere  is  placed  in  the  mixed  layer,  with  sphere  center  at 
15-m  depth,  along  with  the  source  and  receiver  at  10-m 
depth  as  shown  in  Fig.  4(c).  This  is  especially  so  for  the  silt 
bottom,  where  the  scattered  field  from  the  target  sphere  be¬ 
comes  so  greatly  reduced,  when  compared  to  the  previous 
examples  of  this  section,  that  its  returns  only  stand  above  the 
expected  reverberation  level  within  roughly  16-km  range.  It 
is  interesting  that  for  the  sand  bottom,  the  placement  of  the 
target  and  source-receiver  in  the  mixed  layer  leads  to  more 
favorable  conditions  for  detection,  which  should  be  possible 
beyond  50-km  range,  than  if  only  the  source-receiver  were 
placed  in  the  mixed  layer  and  the  target  was  in  the  middle  of 
the  water  column  as  in  Fig.  4(b).  This  is  because  the  higher- 
order  modes  stimulated  by  the  shallow  source,  receiver,  and 
target  can  be  supported  by  the  high-critical-angle  sand  bot¬ 
tom. 

The  exercise  of  changing  the  depths  of  the  source- 
receiver  and  target  is  repeated  in  Fig.  5  for  a  constant  sound 
speed  water  column.  For  the  sand  bottom,  the  level  of  the 
field  scattered  from  the  sphere  is  not  affected  significantly  by 
moving  the  source-receiver  and  target  depths.  For  the  silt 
bottom,  however,  a  significant  decrease  in  the  sphere’s  echo- 
return  level  is  found  for  shallow  source-receiver  and  target 
placements.  Apparently,  these  shallow  placements  stimulate 
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FIG.  9.  Same  as  Fig.  8  except  two-layer  bottom  scenario  of  Fig.  6(d)  is  investigated  for  coherent  scattering  from  subseafloor  riverbank.  (a)  Sediment  is 
comprised  of  light  silt  layer  of  hl=  1-m  thickness  over  a  silt  layer  of  h2  =  1-m  thickness  over  a  sand  half-space.  Coherent  riverbank  scattering  is  from  the 
double  interface  of  light  silt  to  silt  to  sand,  (b)  Sediment  is  comprised  of  light  silt  layer  of  h2  =  1-m  thickness  over  a  sand  layer  of  h2=  1-m  thickness  over 
a  silt  half-space.  Coherent  riverbank  scattering  is  from  the  double  interface  of  light  silt  to  sand  to  silt. 


higher-order  modes  that  are  not  supported  by  the  silt  bottom. 
These  results  should  be  compared  to  those  found  in  Fig.  3(a) 
for  source-receiver  and  target  depths  in  the  middle  of  the 
waveguide  at  50  m. 

B.  Geological  clutter  versus  diffuse  seafloor 
reverberation 

Geomorphic  features  of  the  seafloor  can  return  echoes 
that  stand  well  above  the  diffuse  reverberation  background 
described  in  the  previous  section.  Since  these  echoes  appear 
as  discrete  events  in  time  or  range,  they  may  be  used  to 
remotely  image  seafloor  or  subseafloor  geomorphology  in 
geophysical  applications.  They  may,  however,  also  be  con¬ 
fused  with  returns  from  a  submerged  target  in  an  active  de¬ 
tection  scenario. 

Both  coherent  and  incoherent  scattering  from  the  ca¬ 
nonical  seafloor  and  subseafloor  features,  shown  in  Fig.  6, 
are  investigated.  Both  kinds  of  features  are  modeled  as  a  flat 
100Xl00-m2  surfaces  at  an  inclination  of  10  degrees.  The 
dimensions  and  inclination  are  based  on  actual  geophysical 
data  characterizing  seafloor  and  subseafloor  riverbanks.32 
Seafloor  and  subseafloor  river  channels  are  commonly  found 
in  continental  shelf  waters  after  a  sea  level  rise.  The  latter 
requires  an  additional  influx  of  sedimentation.  In  all  cases  to 
be  considered  here,  the  waveguide  is  modeled  as  an  iso¬ 
velocity  water  column  overlying  one  or  two  sediment  layers 
that  cover  a  sediment  half  space. 

For  the  coherent  calculation,  the  riverbank  is  treated  as  a 
smooth  but  finite  square  surface  with  reflection  coefficient 
appropriate  to  the  given  boundary  conditions,  including  mul¬ 
tiple  reflection  from  various  layers.  The  coherent  scattered 
field  from  the  riverbank  in  the  layered  waveguide  follows 
when  the  scatter  function  for  the  smooth  riverbank,  given  in 
Eq.  (59),  is  inserted  into  Eq.  (1).  Coherent  scattering  from 
the  riverbank  is  then  completely  determined  by  the  boundary 
conditions  at  the  riverbank  and  the  riverbank  geometry.  For 
the  incoherent  calculation,  the  riverbank  is  modeled  first  as  a 


diffusely  scattering  Lambertian  surface  with  the  empirically 
derived  Mackenzie  albedo  and  riverbank  tilt  angle  incorpo¬ 
rated  as  indicated  in  Eq.  (41).  The  diffuse  calculations  are 
also  made  using  perturbation  theory  by  substituting  Eq.  (58) 
into  Eq.  (17).  In  both  cases,  the  assumption  is  that  the  sea¬ 
floor  feature  falls  within  the  resolution  footprint  of  the  sonar 
system. 

Illustrative  examples  are  given  in  Figs.  7-9,  13,  19,  and 
20.  The  geometry  is  again  monostatic  with  colocated  omni¬ 
directional  point  source  and  receiving  array  centers  at  50-m 
depth.  The  receiving  array  lies  parallel  to  the  y  axis.  The 
square  riverbank  surface  has  two  edges  parallel  to  the  y  axis, 
is  centered  at  y  =  0,  and  inclined  10  degrees  about  the  x  axis. 
All  plots  give  the  scattered  field  from  the  riverbank  as  a 
function  of  range  from  the  monostatic  sonar.  For  compari¬ 
son,  incoherent  reverberation  from  the  water-sediment  inter¬ 
face  within  the  resolution  footprint  of  the  sonar,  based  on  the 
Lambert-Mackenzie  model  for  an  un-inclined  surface,  is 
also  plotted  as  a  function  of  range  in  Figs.  7-9  and  13  and 
based  on  perturbation  theory  in  Figs.  19  and  20.  This  is 
referred  to  as  diffuse  seafloor  reverberation.  The  range  and 
cross-range  resolution  of  the  sonar  system  resolution  foot¬ 
print  are  the  same  as  those  stated  in  the  introduction  to  Sec. 
IV.  Similarly,  incoherent  reverberation  from  the  sediment 
layer  to  sediment  half-space  interface,  based  upon  the 
Lambert-Mackenzie  model  for  an  uninclined  surface,  is  also 
plotted  as  a  function  of  range  in  Figs.  7-9  and  13,  and  based 
on  perturbation  theory  in  Figs.  19  and  20.  This  is  referred  to 
as  diffuse  subseafloor  reverberation.  The  far  field  of  the  co¬ 
herent  riverbank  begins  at  roughly  2  km  while  the  far  field  of 
the  receiving  array  begins  at  roughly  1.2  km. 

For  the  Pekeris  waveguide  scenario  of  Fig.  6(a),  returns 
from  the  seafloor  riverbank  features  stand  well  above  diffuse 
seafloor  reverberation  from  the  silt  bottom  within  ranges  of 
roughly  20  km  and  from  the  sand  bottom  beyond  ranges  of 
50  km  when  the  riverbank  is  treated  as  a  coherent  scatterer, 
as  shown  in  Fig.  7.  The  ordinate  is  in  decibels,  i.e.,  20  log|(l)J 
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FIG.  10.  Magnitudes  of  the  coherent  scattering  functions  \S(a,/3=  tt, a t ,  /3,  =  0,)|,  i.e.,  20  log  |S|  dB,  for  the  100X100-m2  seafloor  and  subseafloor  riverbank 
features  at  inclination  \=  10  degrees  of  Fig.  6  over  bistatic  horizontal  grazing  angle  n/2—  a,  for  the  incident  and  a—  7?/2  for  the  scattered  wave,  as  appropriate 
for  backscatter  in  a  waveguide.  The  boxes  include  all  modes  n  where  0.5  rad/km  >Im{f„}.  This  includes  all  and  only  trapped  modes  for  the  Pekeris  waveguide 
scenario  of  Fig.  6(a).  (a)  Reflection  coefficient  for  water  to  sand  is  used  for  scenario  of  Fig.  6(a).  (b)  Reflection  coefficient  of  silt  to  sand  is  used  for  scenario 
of  Fig.  6(c).  (c)  Double  reflection  coefficient  of  water  to  2-m  silt  layer  over  sand  is  used  for  scenario  of  Fig.  6(b).  (d)  Double  reflection  coefficient  of  water 
to  5-m  silt  layer  over  sand  is  used  for  scenario  of  Fig.  6(b).  (e)  Double  reflection  coefficient  of  light  silt  to  1-m  silt  layer  over  sand  is  used  for  scenario  of  Fig. 
6(d). 


for  riverbank  returns  and  10  log  V B  for  diffuse  reverberation. 

For  the  single-layered  bottom  scenarios  of  Figs.  6(b)  and 
(c),  returns  from  both  the  seafloor  riverbank  and  subseafloor 
riverbank  features  can  stand  well  above  diffuse  seafloor  re¬ 


verberation  from  the  silt-over-sand  bottom  when  the  river¬ 
bank  is  treated  as  a  coherent  scatterer  and  the  silt  layer  is  2  m 
or  |  of  a  wavelength,  as  shown  in  Figs.  8(a)  and  (b).  The 
subseafloor  riverbank  and  seafloor  riverbank  return  echoes  at 
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FIG.  1 1 .  Same  as  Fig.  10  except  (a)  reflection  coefficient  for  water  to  silt  is  used  for  scenario  of  Fig.  6(a).  (b)  Reflection  coefficient  of  sand  to  silt  is  used  for 
scenario  of  Fig.  6(c).  (c)  Double  reflection  coefficient  of  water  to  2-m  sand  layer  over  silt  is  used  for  scenario  of  Fig.  6(b).  (d)  Double  reflection  coefficient 
of  water  to  5-m  sand  layer  over  silt  is  used  for  scenario  of  Fig.  6(b).  (e)  Double  reflection  coefficient  of  light  silt  to  1-m  sand  layer  over  silt  is  used  for  scenario 
of  Fig.  6(d). 


similar  levels  within  roughly  5  km,  where  both  features  typi¬ 
cally  stand  above  the  diffuse  seafloor  reverberation  by 
roughly  10  dB,  which  exceeds  the  5.6-dB  standard  deviation. 
The  prominence  of  the  subseafloor  riverbank  returns  follows 
from  the  greater  impedance  mismatch  between  the  silt-sand 
interface  than  the  water-silt  interface  incorporated  in  the  ri¬ 
verbank  scatter  function.  Beyond  roughly  5  km,  the  subsea¬ 


floor  feature  has  returns  that  fall  off  more  rapidly  than  those 
of  the  seafloor  feature.  This  follows  from  the  stripping  of 
higher-order  modes  that  propagate  with  high  attenuation  in 
the  silt  layer.  Coherent  returns  from  the  riverbank  arise  be¬ 
cause  of  its  finite  extent.  Since  the  riverbank  is  modeled  as  a 
smooth  flat  surface,  scattering  is  greatest  in  the  specular  di¬ 
rection  and  falls  off  in  other  directions  in  a  manner  similar  to 
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FIG.  12.  The  horizontal  wave  number’s  imaginary  component  Im{f„}  is  plotted  as  a  function  of  horizontal  grazing  angle,  |ir/2—  a,-|,  for  the  various 
waveguides  considered.  Proper  modes  occur  in  Pekeris  below  the  critical  angle  for  0.5  rad/km  >Im{— £„}.  (a)  Pekeris  with  sand  bottom  and  Pekeris  with  silt 
bottom,  (b)  Constant  water  column  sound  speed  of  1500  m/s  over  It  i  =  2 - m  and  /i[  =  5-m  silt  layer  over  sand  half-space,  (c)  Constant  water  column  sound 
speed  of  1 500  m/s  over  h  j  =  2  -  m  and  h  j  =  5  -  m  sand  layer  over  silt  half-space,  (d)  Constant  water  column  sound  speed  of  1 500  m/s  over  1  -m  light  silt  layer 
over  1-m  silt  layer  over  sand  half-space  and  constant  water  column  sound  speed  of  1500  m/s  over  1-m  light  silt  layer  over  1-m  sand  layer  over  silt  half-space. 


the  sidelobes  of  a  phased  array’s  beampattern.  In  backscatter, 
for  the  given  geometry,  the  riverbank  returns  increase  in  in¬ 
tensity  with  the  square  of  its  length,  or  cross-range  extent. 
Longer  riverbanks  that  fit  within  the  sonar  resolution  foot¬ 
print  then  yield  significantly  larger  returns  as  a  consequence 
of  the  coherent  scattering  assumption,  and  may  stand  well 
above  diffuse  seafloor  reverberation  beyond  10  km.  Returns 
from  such  extended  riverbanks,  however,  rapidly  become 
more  of  a  challenge  to  model  since  the  near  field  moves  out 
in  range  from  the  feature  with  the  square  of  its  length.  As  the 
range  extent  of  the  riverbank  increases,  the  coherent  area 
increases  but  the  side  lobe  level  decreases  for  the  present 
geometry,  rendering  the  effect  on  the  backscattered  field  less 
apparent  than  in  cross-range  augmentation. 

When  the  thickness  of  the  silt  layer  is  increased  to  5  m, 
or  one  wavelength,  returns  from  the  subseafloor  riverbank 
features  are  somewhat  reduced,  as  shown  in  Fig.  8(c),  and 
again  only  stand  above  diffuse  seafloor  reverberation  within 
roughly  5  km.  This  follows  from  a  related  increase  in  the 
stripping  of  the  higher-order  modes  that  have  propagating 
components  in  the  sediment  layer  since  the  sediment  layer 


has  much  higher  attenuation  than  the  water  column.  The  sea¬ 
floor  riverbank  feature  stands  above  diffuse  seafloor  rever¬ 
beration  beyond  20  km  but  rarely  in  excess  of  the  Gaussian 
field  standard  deviation  of  5.6  dB. 

When  the  sediment  layer  is  composed  of  sand  and  the 
half-space  below  is  made  of  silt,  the  situation  changes  dras¬ 
tically,  as  shown  in  Figs.  8(d)  and  (e).  Returns  from  the 
subseafloor  riverbank  no  longer  stand  above  diffuse  seafloor 
reverberation  beyond  2  km  because  in  the  sand  layer,  which 
is  much  faster  than  the  water  column  and  silt  half-space,  the 
trapped  modes  become  evanescent.  Seafloor  riverbank  re¬ 
turns  stand  well  above  diffuse  seafloor  reverberation,  occa¬ 
sionally  by  10  dB  or  more,  even  beyond  20  km  for  both  the 
2-  and  5-m-thick  sand  layers,  as  expected  given  the  large 
impedance  contrast  between  water  and  sand.  Older  sea¬ 
floor  features,  in  fact,  are  more  likely  to  be  composed  of 
consolidated  material  such  as  sand  or  limestone  since  such 
materials  are  better  able  to  withstand  erosion.  Steeper  sea¬ 
floor  features  that  are  common  in  many  continental  shelves, 
such  as  glacier  and  iceberg  scours,  can  yield  even  higher 
returns. 
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FIG.  13.  Same  as  Figs.  7  and  8  except  the  Lambert-Mackenzie  model  is  used  to  model  scattering  from  inclined  riverbank  features,  (a)  Seafloor  riverbank  over 
sand  and  silt  half-spaces,  (b)  Seafloor  and  subseafloor  riverbank  scattering  with  the  upper  sediment  layer  composed  of  silt  with  h  x  =  2  m  and  the  lower 
sediment  half-space  composed  of  sand,  (c)  Same  as  (b)  except  hl  =  5  m.  (d)  Seafloor  and  subseafloor  riverbank  scattering  with  the  upper  sediment  layer 
composed  of  sand  with  h  j  =  2  m  and  the  lower  sediment  half-space  composed  of  silt,  (e)  Same  as  (d)  except  h  { =  5  m. 


For  the  two-layered  bottom  of  Fig.  6(d),  the  subseafloor 
riverbank  returns  stand  roughly  10  dB  above  diffuse  seafloor 
reverberation  out  to  roughly  10  km  for  a  1-m  light-silt  layer 
over  a  1-m  silt  layer  over  a  sand  half-space,  as  shown  in  Fig. 
9(a).  The  double  layer  reflection  coefficient  from  the  light- 
silt  to  silt  to  sand  interfaces  leads  to  the  increased  promi¬ 
nence  of  the  subseafloor  riverbank  returns,  compared  with 
those  obtained  with  the  single -layer  reflection  coefficient  of 


Figs.  8(a)-(e).  When  the  layering  is  altered  to  1-m  light-silt 
over  sand  over  a  silt  half-space,  returns  from  the  subseafloor 
riverbank  feature  only  stand  above  diffuse  seafloor  rever¬ 
beration  within  roughly  5  km.  This  indicates  that  sediment 
stratification  of  the  geomorphic  feature  can  weigh  in  heavily 
in  fixing  its  scattering  amplitude. 

The  effect  of  bottom  layering  on  the  coherent  scattering 
function  of  the  inclined  seafloor  and  subseafloor  riverbank 
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FIG.  14.  Scattering  strength  SS(0=  77— 0, ,  </>=0;  0,  ,  </>,=  77)  in  free-space 
backscatter  as  a  function  of  surface  grazing  angle  |  it  12—  0,|  for  a  diffusely 
scattering  surface  obeying  Lambert-Mackenzie  and  first-order  perturbation 
theory  scattering  laws.  The  first-order  perturbation  theory  curves  are  for 
cases  where  the  plane  wave  is  incident  from  an  upper  to  a  lower  medium, 
where  the  upper  and  lower  media  can  be  water,  sand,  or  silt. 

features  of  Fig.  6  is  illustrated  in  Figs.  10  and  11  as  a  func¬ 
tion  of  horizontal  grazing  angle,  7r/2—  a,  for  the  incident  and 
a—  77-/2  for  the  scattered  wave,  at  fixed  incident  and  scattered 
azimuths,  /3,  =  0,  (3=tt,  as  is  appropriate  to  backscatter  in  a 
waveguide.  The  trapped  modes  for  the  Pekeris  waveguide 
scenarios,  Fig.  10(a)  for  a  sand  bottom  and  Fig.  11(a)  for  a 
silt  bottom,  have  incident  and  scattered  elevation  angles  that 
lie  within  the  boxes  shown,  the  dimensions  of  which  corre¬ 
spond  to  the  respective  bottom  critical  angle.  The  boxes  in¬ 
clude  all  modes  n  where  0.5  rad/km>Im{£„}.  This  includes 
all  and  only  trapped  modes  for  the  Pekeris  waveguide  sce¬ 
nario  of  Fig.  6(a).  The  latter  criterion  is  used  to  segment 
modes  that  dominate  the  incident  propagation  by  similar 
boxes  for  the  more  complicated  layered  bottom  cases  illus¬ 
trated  in  Figs.  10  and  11.  The  value  Im{£„}  is  plotted  as  a 
function  of  equivalent  modal  angle  a  in  Fig.  12  for  the  vari¬ 
ous  waveguides  considered.  This  makes  it  possible  to  see 
how  the  scattering  functions  of  Figs.  10  and  11  are  discretely 
sampled  in  the  waveguide  scattering  theory  defined  by  Eq. 
(1)  and  to  estimate  the  attenuation  of  a  given  modal  compo¬ 
nent  as  a  function  of  range.  Inspection  of  Figs.  10  and  11 
reveals  that  seafloor  and  subseafloor  riverbank  features  that 
backscatter  most  prominently  in  Figs.  7-9  have  scatter  func¬ 
tions  with  relatively  large  amplitudes  at  the  equivalent  angles 
of  the  propagating  modes.  While  modes  propagating  at 
steeper  angles  suffer  greater  attenuation,  as  indicated  in  Fig. 
12,  these  same  modes  are  scattered  much  more  efficiently  by 
the  slightly  inclined  riverbank  features  as  indicated  in  Figs. 
10-11,  so  that  there  is  some  balancing  between  the  two  ef¬ 
fects  that  is  unique  to  waveguide  scattering.  Higher-order 
modes,  with  elevation  angles  less  than  45  degrees  where  0 
degrees  points  downward,  however,  contribute  negligibly  to 
the  field  scattered  from  the  riverbank  features  for  the  ranges 
and  features  investigated  in  the  present  article. 

When  the  riverbank  is  treated  as  an  incoherent  scatterer 
with  the  Lambert-Mackenzie  model  of  Eq.  (41),  returns 


from  both  the  100X100-m2  seafloor  and  subseafloor  river¬ 
bank  features  at  10  degrees  inclination  never  stand  above 
diffuse  seafloor  reverberation  by  more  than  a  fraction  of  the 
expected  5.6-dB  standard  deviation,  as  shown  in  Fig.  13. 
Riverbank  returns  are  again  in  decibels,  i.e.,  10  times  the  log 
of  the  covariance  given  in  Eq.  (41).  This  is  still  the  case  for 
the  ranges  shown  in  Fig.  13,  except  for  the  2-m  sand  layer, 
even  if  the  riverbank  feature  is  extended  laterally  to  fill  the 
entire  cross-range  width  of  the  system  resolution  footprint,  as 
can  be  readily  checked  by  noting  that  diffuse  reverberation 
accrues  in  direct  proportion  to  the  area  of  the  scattering 
patch.  For  the  2-m  sand  layer,  the  seafloor  riverbank  can 
have  returns  that  exceed  the  diffuse  reverberation  back¬ 
ground  by  more  than  5.6  dB  if  it  fills  the  entire  resolution 
footprint  in  cross-range. 

Diffuse  subbottom  reverberation,  shown  in  Fig.  13,  al¬ 
ways  returns  at  a  lower  level  than  diffuse  seafloor  reverbera¬ 
tion  if  the  same  empirical  Lambert-Mackenzie  incoherent 
scattering  law  is  used.  This  comparison  highlights  the  differ¬ 
ences  in  propagation  to  and  from  the  seafloor  and  subbottom 
interfaces  since  the  scattering  function  is  held  fixed.  The 
comparison  may  be  purely  academic,  however,  because  the 
Lambert-Mackenzie  law  serves  as  an  empirical  catch-all  that 
describes  the  entire  seabed  scattering  process  and  so  already 
incorporates  the  effect  of  bottom  layering  and  volume  scat¬ 
tering  in  some  average  sense.  There  is,  in  other  words,  no 
reason  to  believe  that  scattering  from  the  different  interfaces 
can  be  modeled  with  exactly  the  same  albedo  and  scattering 
law. 

Perturbation  theory  offers  a  more  fundamental  approach 
to  modeling  rough  surface  scattering  that  can  also  be  used  to 
investigate  potential  mechanisms  for  geological  clutter. 
While  the  impedance  contrast  at  the  scattering  interface  is 
fully  accounted  for  in  the  perturbation  theory  formulation, 
additional  parameters  describing  the  roughness  spectrum 
must  be  known.  The  perturbation  theory  formulation  de¬ 
scribed  in  Sec.  IIID  is  used  with  the  spectral  strength  and 
power  law  parameters  W2  =  0. 04/(2  7r)  and  y=4.0,  yielding 
frequency-independent  scattering,  following  Essen.  ’3  These 
values  are  not  based  on  physical  measurements,  since  none 
are  presently  available  in  the  present  frequency  range,  but 
rather  have  been  chosen  so  that  the  scattering  strength  that 
perturbation  theory  yields  is  near  that  of  the  empirical 
Lambert-Mackenzie  model  for  the  various  single  and  mul¬ 
tiple  reflection  interfaces  considered  here,  as  shown  in  Fig. 
14.  In  all  curves  where  scattering  arises  from  a  wave  incident 
from  a  slower  medium,  a  discontinuity  in  slope  is  found  at 
the  critical  angle.  Beyond  this  a  significant  reduction  in  scat¬ 
tering  occurs  until  roughly  45  degrees  where  shallow  angle 
assumptions  of  first-order  perturbation  theory  are  no  longer 
valid  for  the  given  surface  roughness  parameters,  and  the 
curves  increase  dramatically  in  an  unphysical  manner.  Dif¬ 
ferences  in  the  perturbation  theory  curves  away  from  the 
critical  angles  arise  principally  from  the  impedance  contrasts 
between  the  media  considered.  An  exception  occurs  for  fast 
sand  over  slow  silt  where  no  critical  angle  exists  and  trans¬ 
mission  into  the  silt  is  significant  even  at  very  shallow 
angles,  where  low  level  scattering  results. 

The  effect  of  scattering  and  reflection  from  multiple  lay- 
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FIG.  15.  Scattering  strength  SS  1  cr.  / 3=tt ,  a , ,  /3,  =  0,),  based  on  first-order  perturbation  theory  for  the  seafloor  and  subseafloor  riverbank  features  at  inclination 
X=  10  degrees  of  Fig.  6  over  bistatic  horizontal  grazing  angle  77/2—  or,  for  incident  waves  and  a—  tt/2  for  scattered  waves,  as  appropriate  to  backscatter  in  a 
waveguide.  The  boxes  include  all  modes  n  where  0.5  rad/km  >Im{f„}.  This  includes  all  and  only  trapped  modes  for  the  Pekeris  waveguide  scenario  of  Fig. 
6(a).  (a)  Reflection  coefficient  for  water  to  sand  is  used  for  scenario  of  Fig.  6(a).  (b)  Reflection  coefficient  of  silt  to  sand  is  used  for  scenario  of  Fig.  6(c).  (c) 
Double  reflection  coefficient  of  water  to  2-m  silt  layer  over  sand  is  used  for  scenario  of  Fig.  6(b).  (d)  Double  reflection  coefficient  of  water  to  5-m  silt  layer 
over  sand  is  used  for  scenario  of  Fig.  6(b).  (e)  Double  reflection  coefficient  of  light  slit  to  1-m  silt  layer  over  sand  is  used  for  scenario  of  Fig.  6(d). 


ers  can  be  significant  as  shown  in  Figs.  15  and  16  for  an 
inclined  riverbank  surface  and  Figs.  17  and  18  for  general 
uninclined  seafloor,  where  the  perturbation  theory  scattering 
strength  is  presented  for  the  bistatic  scattering  scenario  rel¬ 


evant  to  backscatter  in  a  waveguide,  as  was  done  in  Figs.  10 
and  1 1  for  the  coherent  scatter  function. 

Empirical  values  for  the  spectral  strength  and  power  law 
parameters  of  first-order  perturbation  theory  have  been  ob- 
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FIG.  16.  Same  as  Fig.  15  except  (a)  reflection  coefficient  for  water  to  silt  is  used  for  scenario  of  Fig.  6(a).  (b)  Reflection  coefficient  of  sand  to  silt  is  used  for 
scenario  of  Fig.  6(c).  (c)  Double  reflection  coefficient  of  water  to  2-m  sand  layer  over  silt  is  used  for  scenario  of  Fig.  6(b).  (d)  Double  reflection  coefficient 
of  water  to  5-m  sand  layer  over  silt  is  used  for  scenario  of  Fig.  6(b).  (e)  Double  reflection  coefficient  of  light  silt  to  1-m  sand  layer  over  silt  is  used  for  scenario 
of  Fig.  6(d). 


tained  for  various  seafloor  types  by  Jackson34  over  the  short 
spatial  scales  relevant  to  the  analysis  of  high-frequency  scat¬ 
tering  in  the  ten  kilohertz  range  and  beyond,  where  this  au¬ 
thor  has  shown  perturbation  theory  to  match  experimental 
data  well.  When  these  same  values  are  used  at  low  fre¬ 
quency,  specifically  /=  300  Hz,  the  resulting  scattering  law 
falls  more  than  an  order  of  magnitude  below  the  empirical 


Lambert-Mackenzie  curve  shown  in  Fig.  14.  Unrealistically 
high  roughness  values  for  the  spectral  strength,  as  obtained 
for  rough,  rocky  surfaces  in  the  high-frequency  analysis  of 
Ref.  34,  are  necessary  for  first-order  perturbation  theory  to 
match  the  empirical  seafloor  scattering  strength  curve  of 
Mackenzie  in  the  low-frequency  regime  of  interest  here. 
Since  the  Mackenzie  curve  summarizes  the  entire  seafloor 
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FIG.  17.  Scattering  strength  SS(a,  P=ir,  a , ,  /3,  =  0,)  based  on  first-order  perturbation  theory  for  level  seafloor,  ^=0  degrees,  over  bistatic  horizontal  grazing 
angle  tt/2—  ai  for  the  incident  wave  and  a—  ir/2  for  the  scattered  wave,  as  appropriate  to  backscatter  in  a  waveguide.  The  boxes  include  all  modes  n  where 
0.5>Im{f„}.  This  includes  all  and  only  trapped  modes  for  the  Pekeris  waveguide  scenario  of  Fig.  6(a).  (a)  Reflection  coefficient  for  water  over  sand  is  used, 
(b)  Double  reflection  coefficient  of  water  over  2-m  silt  over  sand  is  used,  (c)  Double  reflection  coefficient  of  water  over  5-m  silt  layer  over  sand  is  used,  (d) 
Triple  reflection  coefficient  of  water  over  1-m  light  silt  layer  over  1-m  silt  layer  over  sand  is  used. 


scattering  process,  and  is  not  limited  to  interface  scattering 
or  the  more  restrictive  type  of  interface  scattering  described 
by  first-order  perturbation  theory,  it  is  reasonable  to  conclude 
that  either  the  assumptions  of  first-order  perturbation  theory 
are  inadequate  to  properly  model  seafloor  scattering  at  low 
frequency,  a  significantly  different  set  of  spectral  strength 
and  power  law  parameters  must  characterize  seafloor  inter¬ 
face  scattering  at  low  frequency,  or  a  more  sophisticated 
modeling  of  the  seabed  layering  and  sound  speed  gradients  is 
necessary.  It  is  also  possible  that  scattering  from  volume 
heterogeneities  may  yield  significant  reverberation.  This  is 
most  likely  to  be  the  case  where  a  propagating,  rather  than 
evanescent,  component  of  the  modal  spectrum  exists  in  the 
layer  where  the  volume  heterogeneities  are  present. 

When  the  riverbank  feature  is  treated  as  an  incoherent 


scatterer  using  first-order  perturbation  theory,  as  in  Figs.  19 
and  20,  only  seafloor  riverbank  returns  from  a  single  sand 
layer  can  stand  above  diffuse  seafloor  reverberation  by  more 
than  5.6  dB,  and  this  only  occurs  when  the  riverbank  feature 
is  extended  laterally  to  fill  the  entire  cross-range  extent  of  the 
system  resolution  footprint.  Returns  from  the  subseafloor  ri¬ 
verbank  only  stand  above  diffuse  seafloor  reverberation  by 
more  than  5.6  dB  for  the  two-layered  bottom  in  the  light-silt 
over  sand  over  silt  scenario,  and  this  only  occurs  if  the  fea¬ 
ture  is  extended  to  fill  the  resolution  footprint. 

Comparison  of  Figs.  7  and  9  and  Fig.  13  shows  that 
coherent  returns  greatly  outweigh  incoherent  returns  from 
the  riverbank  feature.  This  finding  is  advantageous  since 
only  deterministic  physical  and  geometrical  parameters  of 
the  seafloor  are  necessary  in  the  coherent  model,  whereas 
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FIG.  18.  Same  as  Fig.  17  except  (a)  reflection  coefficient  for  water  over  silt  is  used,  (b)  Double  reflection  coefficient  of  water  over  2-m  sand  over  silt  is  used, 
(c)  Double  reflection  coefficient  of  water  over  5-m  sand  layer  over  silt  is  used,  (d)  Triple  reflection  coefficient  of  water  over  1-m  light  silt  layer  over  1-m  sand 
layer  over  silt  is  used. 


either  empirical  data  or  a  stochastic  representation  of  the 
seafloor  is  necessary  in  the  incoherent  model.  The  environ¬ 
mental  description  necessary  for  the  coherent  model  is  then 
easier  to  obtain  and  rests  on  far  fewer  supporting  assump¬ 
tions  than  the  incoherent  one. 

V.  CONCLUSIONS 

One  of  the  greatest  challenges  to  active  sonar  operations 
in  shallow  water  arises  when  echo  returns  from  the  intended 
target  become  indistinguishable  from  reverberation  returned 
by  the  waveguide  boundaries  and  volume.  To  determine  con¬ 
ditions  in  which  a  typical  low-frequency  active  sonar  system 
may  operate  effectively  in  a  shallow  water  waveguide,  a  uni¬ 
fied  model  for  submerged  object  scattering  and  reverberation 
is  developed.  The  approach  is  to  use  a  waveguide  scattering 
model  that  follows  directly  from  Green’s  theorem  but  that 
takes  advantage  of  simplifying  single-scatter  and  far-field  ap¬ 
proximations  that  apply  to  a  wide  variety  of  problems  where 


the  source  and  receiver  are  distant  from  the  target.  To  treat 
reverberation  from  randomly  rough  boundaries  and  stochas¬ 
tic  volume  inheterogeneties,  the  waveguide  scattering  model 
is  generalized  to  include  stochastic  targets.  Analytic  expres¬ 
sions  for  the  spatial  covariance  of  the  field  scattered  from  a 
stochastic  target  are  then  obtained  in  terms  of  the  waveguide 
Green’s  function  and  the  covariance  of  the  target’s  plane 
wave  scatter  function.  This  makes  the  formulation  amenable 
to  a  wide  variety  of  approaches  for  computing  a  stochastic 
target’s  scatter  function.  For  diffuse  seafloor  reverberation, 
two  approaches  are  adopted,  an  empirical  one  of  Lambert 
and  Mackenzie  and  a  fundamental  one  based  on  first-order 
perturbation  theory.  It  is  most  convenient  to  describe  the 
diffuse  component  of  distant  seafloor  reverberation  with  a 
modal  formulation  since  the  modes  comprise  the  statistical 
entities  of  the  field  that  the  scattering  surface  may  decorre¬ 
late. 

Since  reverberation  is  measured  in  time  but  the  wave¬ 
guide  scattering  formulation  is  for  harmonic  field  compo- 
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FIG.  19.  Same  as  Figs.  7  and  8  except  first-order  perturbation  theory  is  used  to  model  scattering  from  the  inclined  riverbank  features,  (a)  Seafloor  riverbank 
over  sand  and  silt  half-spaces,  (b)  Seafloor  and  subseafloor  riverbank  scattering  with  the  upper  sediment  layer  composed  of  silt  with  hl  =  2  m  and  the  lower 
sediment  half-space  composed  of  sand,  (c)  Same  as  (b)  except  hl  =  5  m.  (d)  Seafloor  and  subseafloor  riverbank  scattering  with  the  upper  sediment  layer 
composed  of  sand  with  h1=2  m  and  the  lower  sediment  half-space  composed  of  silt,  (e)  Same  as  (d)  except  hl  =  5  m. 


nents,  the  time  dependence  of  the  field  scattered  by  a  distant 
object  from  a  source  of  arbitrary  time  dependence  is  derived 
analytically  using  the  saddle  point  method.  The  resulting  ex¬ 
pression  is  given  in  terms  of  modal  group  velocities,  the 
frequencies  of  which  vary  as  a  function  of  time  and  source, 
receiver,  and  target  position.  A  simpler  analytic  approach 
involving  Parseval’s  theorem  can  be  applied  when  the  inte¬ 
gration  time  of  the  measurement  system  is  sufficiently  long 


to  include  the  dominant  energy  returned  from  the  target  or 
scattering  patch.  This  approach  is  used  in  the  illustrative  ex¬ 
amples.  A  viewer-oriented  reference  frame  is  then  adopted, 
translating  from  the  traditional  target-oriented  frame  of 
waveguide  scatter  theory,  to  incorporate  the  continuous  dis¬ 
tribution  of  scatterers  encountered  in  waveguide  boundary 
and  volume  reverberation.  This  enables  analytic  expressions 
to  be  developed  for  the  reverberant  field  returned  bistatically 
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FIG.  20.  Same  as  Fig.  9  except  first-order  perturbation  theory  is  used  to  model  scattering  from  the  inclined  subseafloor  riverbank  for  the  two-layer  bottom 
scenario  of  Fig.  6(d).  (a)  Sediment  is  comprised  of  light  silt  layer  of  h  j  =  1  -m  thickness  over  a  silt  layer  of  h  ,  =  1  -m  thickness  over  a  sand  half-space.  Diffuse 
riverbank  scattering  is  from  the  double  interface  of  light  silt  to  silt  to  sand,  (b)  Sediment  is  comprised  of  light  silt  layer  of  h  l  =  1  -  m  thickness  over  a  sand  layer 
of  li2  =  1-m  thickness  over  a  silt  half-space.  Diffuse  riverbank  scattering  is  from  the  double  interface  of  light  silt  to  sand  to  silt. 


from  seafloor  within  the  resolution  footprint  of  a  typical  ac¬ 
tive  sonar  system  after  narrow-band  beamforming  with  a 
horizontal  array. 

The  unified  model  is  used  to  investigate  typical  low- 
frequency  active  detection  scenarios  in  shallow  water. 
Sample  calculations  for  finite-duration  cw  source  signals  in¬ 
dicate  that  the  maximum  range  at  which  echo  returns  from  a 
submerged  target  stand  unambiguously  above  diffuse  sea¬ 
floor  reverberation  is  highly  dependent  upon  the  water  col¬ 
umn  and  sediment  stratification,  as  well  as  the  receiving  ar¬ 
ray  aperture,  source,  receiver,  and  target  location,  and  the 
scattering  properties  of  the  target  and  seafloor. 

The  model  is  also  applied  to  determine  conditions  in 
which  discrete  morphological  features  of  the  seafloor  and 
subseafloor  return  echoes  that  stand  prominently  above  dif¬ 
fuse  seafloor  reverberation.  Simulations  for  finite-duration 
cw  source  signals  indicate  that  typical  seafloor  and  subsea¬ 
floor  riverbank  features,  ubiquitously  found  throughout  con¬ 
tinental  shelf  waters,  can  return  echoes  that  stand  signifi¬ 
cantly  above  the  diffuse  component  of  seafloor  reverberation 
in  the  operational  ranges  of  typical  low-frequency  active  so¬ 
nar  systems.  This  finding  is  significant  since  returns  from 
these  discrete  features  can  be  confused  with  returns  from  an 
intended  submerged  target.  The  relative  prominence  of  this 
kind  of  geological  clutter  is  highly  dependent  on  the  wave¬ 
guide  properties,  measurement  geometry,  and  scattering 
characteristics  of  the  geological  feature  and  surrounding  sea¬ 
floor.  The  finding  that  subseafloor  features  can  cause  signifi¬ 
cant  clutter  is  particularly  troubling  for  active  sonar  opera¬ 
tions  because  it  greatly  increases  the  environmental 
characterization  necessary  to  make  accurate  predictions  of 
the  expected  clutter.  The  coherent  component  of  the  field 
scattered  from  the  riverbank  features  examined,  arising  from 
the  features’  finite  size,  is  found  to  far  outweigh  the  diffuse 
component  arising  from  random  roughness  of  the  features. 
The  methods  and  findings  of  this  article  are  presently  being 
used  to  help  design  a  number  of  field  experiments  to  inves¬ 
tigate  the  physical  mechanisms  that  lead  to  geological  clutter 


in  the  output  of  active  sonar  systems  operated  in  shallow 
water. 
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APPENDIX  A:  THE  PLANE  WAVE  SCATTER 
FUNCTION,  TARGET  STRENGTH,  SURFACE 
SCATTERING  STRENGTH,  AND  GREEN’S  THEOREM 

Standard  parameters  used  to  describe  surface  and  target 
scattering  in  ocean  acoustics  can  be  traced  back  to  Green’s 
theorem  by  using  some  simple  approximations  involving  the 
plane  wave  scatter  function  in  free  space.  To  do  so,  it  must 
be  recalled  that  the  harmonic  field  d>s(r)  scattered  by  an 
object  can  be  expressed  in  terms  of  the  medium  Green  func¬ 
tion  G(r|r,)  and  incident  field  by  the  Helmholtz  - 

Kirckoff  integral  equation  ’5 

f  f  <?G(r|r,) 

®,(r)  =  -J  J  [<t>;(r,)  +  d>.s(r,)] 

A, 

—  G(r|r,)  -  [d>;(r,)  +  <t>s(r,)](fA,,  (Al) 

ont 

a  form  of  Green’s  theorem,  where  G(r|rf)  and  <I>,( r)  each 
satisfy  the  Helmholtz  equation,  driven  by  a  source  at  angular 
frequency  u>  =  2Ttf.  The  area  integral  encloses  the  scatterer 
and  the  surface  normal  points  into  the  enclosed  volume. 
Consider,  first,  the  problem  of  a  plane  wave 

d>,(r  t)  =  eikr>^e>'^\  (A2) 

incident  on  an  object  in  free  space  traveling  in  the  direction 
( 0j ,  d>j)  where,  for  example, 
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FIG.  Bl.  A  comparison  of  scattering  using  the  single  frequency  approximation  versus  the  full  bandwidth  of  the  given  window  function.  Reverberation 
calculated  using  the  single  frequency  approximation  is  indistinguishable  from  that  calculated  with  the  full  bandwidth,  (a)  Same  as  Fig.  3(a)  for  sphere  in 
waveguide  with  reverb  except  only  the  sand  bottom  case  is  shown.  Single  frequency  approximation  is  compared  to  rectangular  window,  (b)  Same  as  (a)  except 
the  Hamming  window  is  used  instead  of  the  rectangular  window,  (c)  Same  as  (a)  except  only  the  silt  bottom  case  is  shown,  (d)  Same  as  (c)  except  the 
Hamming  window  is  used  instead  of  the  rectangular  window. 


r/{  6j ,  <f>i)  =  cos  6j  cos  #f  +  sin  0t  sin  0,  cos(  </>,■—  </>,). 


The  Green  function  is 

1  e‘k\r-r,\  j  eik[r-rtV(0,<p )] 

477  1 1 —  rf |  477  r 


G(r|r,)  = 


(A3) 


(A4) 


where  the  last  approximation  is  for  a  receiver  so  far  from  the 
object  that  r>rt.  By  application  of  Green’s  theorem,  the 
scattered  field  at  this  distant  receiver  then  can  be  written  as 


, ikr 


J  J  [eikr'v<e'  ’A-'  +  cp^r,)]  ~^—e-ikrtv(0,<t>) 


_e-ikrtV(e,4> ) — iy*r, $  ( r  )-\dAt.  (A5) 

on, 

By  the  definition  of  the  plane  wave  scatter  function 
S(  0,  (f>\  Oj ,  </>,),  however,  it  can  also  be  written  as 


_ ikr 


(A6) 


in  an  object-centered  coordinate  system,  which  leads  to  the 
equality 


= - —  f  f  +  (r  )]  —e-ikr,v(6,<l>) 

477  J  J  J  f  J 


—  e~ikrrV(  6,4>) — +  (r MdAt,  (A7) 

Otlf 


which  relates  Eq.  (A5)  directly  to  Green’s  theorem  when  r 
>rt. 

Equation  (A5)  can  be  recast  as  a  sonar  equation  by  tak¬ 
ing  10  log  of  the  squared  magnitude  of  both  sides.  In  terms 
of  the  plane  wave  scatter  function,  the  resulting  target 
strength  of  the  scatterer  is  then 


^lOlog 


SWfcO,,#,)  2 

k 


dB  re  1  m. 


(A8) 
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FIG.  B2.  A  comparison  of  scattering  using  the  single  frequency  approximation  versus  the  full  bandwidth  of  the  given  window  function.  Reverberation 
calculated  using  the  single  frequency  approximation  is  indistinguishable  from  that  calculated  with  the  full  bandwidth,  (a)  Same  as  Fig.  7(a)  for  seafloor 
riverbank  with  reverb  except  only  the  sand  bottom  case  is  shown.  The  single  frequency  approximation  is  compared  to  the  rectangular  window,  (b)  Same  as 
(a)  except  the  Hamming  window  is  used  instead  of  the  rectangular  window,  (c)  Same  as  (a)  except  only  the  silt  bottom  case  is  shown,  (d)  Same  as  (c)  except 
the  Hamming  window  is  used  instead  of  the  rectangular  window. 


The  differential  scattering  cross  section  of  the  target,  defined  /( 9,  </>;  0, ,  (/>,),  a  standard  empirical  descriptor  of  surface 
as  scattering  properties  in  modern  radiometry,  via 


cr=  lim47rr 


<*>,(r)|2  ’ 


(A9) 


= /( 9, 4>\  Oj ,  <t>j)  cos  9  cos  9/AA, 


(All) 


then  becomes 


cr(  0, 4>\  Oj ,  (f>j)  =  4  tt 


k 


(A  10) 


which  in  the  high-frequency  limit  can  be  interpreted  as  the 
projected  area  of  the  target  as  seen  with  foreshortening  from 
the  combined  perspectives  of  the  source  and  receiver. 

If  the  target  is  a  random  patch  of  rough  surface  rather 
than  a  finite  object,  Eqs.  (A5)  and  (A6)  can  still  be  used  so 
long  as  far-held  conditions  hold  and  S(  9,  </>;  9, ,  (/>,)  is  inter¬ 
preted  as  a  stochastic  parameter.  The  scattering  coefficient  of 
the  surface  is  the  expected  scattering  cross-section  per  solid 
angle  (cr/47r).  This  can  be  written  in  terms  of  the  bidirec¬ 
tional  scattering  (reflectance)  distribution  function29 


where  A ,4  is  the  area  of  the  surface  patch.  For  diffuse  scat¬ 
tering  that  is  equal  in  all  directions,  the  bidirectional  scatter¬ 
ing  distribution  function  equals  the  constant  eh r,  where  e  is 
the  albedo,29  which  is  unity  for  a  perfect  Lambertian  surface 
and  roughly  7r/1027  for  a  Lambert-Mackenzien  surface.9 
The  conventional  seafloor  scattering  strength  of  ocean 
acoustics22  is  then 


SS(9,4>;9j,<f>i)=  10  log' 


SWfaO,,#,) 

k 


-10  log  AA, 

(A  12) 


where  the  first  term  on  the  right-hand  side  can  be  interpreted 
as  an  equivalent  target  strength  for  the  random  scattering 
patch. 
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Many  surfaces  scatter  with  both  a  deterministic  as  well 
as  a  fluctuating  component.  It  is  conventional  to  assume  that 
the  fluctuating  component  scatters  a  zero-mean  field.  Diffuse 
scattering  surfaces,  as  described  by  Lambert’s  law  or  first- 
order  perturbation  theory,  scatter  only  zero-mean  fluctuating 
fields. 

The  preceding  analysis  shows  that  the  standard  param¬ 
eters  used  to  describe  surface  and  target  scattering  in  free 
space  can  be  traced  back  to  Green’s  theorem  through  the 
plane  wave  scatter  function. 


APPENDIX  B:  SINGLE  FREQUENCY  APPROXIMATION 
VERSUS  FULL  BANDWIDTH  IN  NARROW-BAND 
SCATTERING  CALCULATIONS 

A  single  frequency  approximation  is  used  for  all  the 
narrow-band  scattering  calculations  of  Sec.  IV.  The  time- 
averaged,  expected  mutual  intensity  of  Eq.  (35), 

7(R,R,  ,R()  ,t)=l-  Jje(/)|2<  |  ®,(Rr  -  R|  Rq  -  R)  1 2)df, 

(Bl) 

is  approximated  as 

/(R,R,.,R0,0~(|<5j(R,-R|Ro-R)|2> 

jjQ{f)\2df,  (B2) 

where  (|<i>s(R(.  —  R|R()  —  R)|2)  is  calculated  at  the  center  fre¬ 
quency  fc  of  the  narrow-band  waveform  q(t),  where  fc 
=  300  Hz  for  the  examples  of  Sec.  IV.  In  practice,  the 
narrow-band  waveform  or  window  function  q(t)  has  a  spec¬ 
trum  Q{f)  with  either  a  narrow  main  lobe  such  as  the  rect¬ 
angular  window,  or  a  broader  main  lobe  with  lower  side 
lobes  such  as  the  Bartlett,  Hanning,  and  Hamming 
windows.36  The  window  functions  are  normalized  according 
to 

\\[\Q(f)\2df=l,  (B3) 


so  that  the  rectangular  window,  with  a  main  lobe  half-power 
bandwidth  of  0.886/7  and  a  first  side  lobe  13.4  dB  down 
peak-to-peak  from  the  main  lobe,  becomes 


q(t) 


1 

0 


for  -  772  772  , 

otherwise. 


(B4) 


and  the  Hamming  window,  for  example  with  main  lobe  half¬ 
power  bandwidth  of  1.30 IT  and  a  first  side  lobe  42.7  dB 
down  peak-to-peak  from  the  main  lobe,  becomes 


q(t)  = 


Vl/0.3974{0.54+0.46  cos(2  irt/T)} 
for  -  7/2  7/2, 

0  otherwise, 


(B5) 


where  T=  1/2  s  for  the  examples  of  Sec.  IV. 

For  the  reverberation  calculations  of  Sec.  IV,  calcula¬ 
tions  show  that  this  approximation  differs  from  the  full  spec¬ 
tral  integration  by  less  than  0.1  dB. 

As  may  be  expected  in  the  coherent  scattering  from  tar¬ 
gets  where  modal  interference  is  significant,  some  range- 


dependent  nulls  and  valleys  in  the  sound  pressure  level  of  the 
received  field  found  in  the  single  frequency  calculation  may 
be  partially  filled  when  the  full  bandwidth  is  used  for  the 
narrow-band  waveforms  considered.  This  is  exhibited  in 
Figs.  Bl  and  B2  where  the  filling  is  shown  to  be  window 
dependent  and  more  negligible  for  bottoms  that  support 
fewer  trapped  modes.  It  is  noteworthy  that  for  narrow-band 
transmissions  at  the  given  center  frequency  and  duration,  the 
sphere  target  may  have  returns  that  fall  below  the  expected 
reverberation  level,  but  the  quantity  and  location  of  these 
expected  deep  “fades”  of  the  target  are  highly  dependent  on 
the  window  function  used. 
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